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ABSTRACT 


An investigation of the interaction effects of a lower heated tube on pool boiling of 
pure R-124 from an upper horizontal tube was conducted at a saturation temperature of 
2.2?C. The test tubes used were: (1) smooth tubes and (2) deformed surface (TURBO-B) 
enhanced tubes. The effects of tube spacing/configuration and lower tube heat flux on the 
heat transfer performance of the upper tube were investigated. 

For both tube arrays, the enhancing effect of bubbles from a lower tube was 
dramatic. This enhancement increased as lower tube heat flux increased. However, when 
upper tube heat fluxes were greater than 20 kW/m^, all enhancement disappeared. 

For a smooth tube array in natural convection, the effect of a lower heated tube on 
the heat transfer from an upper tube was small. In nucleate boiling, a P/D of 1.8 gave the 
best upper tube heat transfer performance and a vigorously nucleating lower tube 
eliminated upper tube hysteresis. With the lower tube unheated and an upper tube heat 
flux of greater than 3 kW/m^, the performance using R-124 was generally better than for 
R-114. With a nucleating lower tube (at 10 kW/m’), again the performance of R-124 was 
better, but only for upper tube heat fluxes of greater than 40 kW/m’. 

For a TURBO-B tube array, a 30 degree offset of the upper tube reduced the upper 
tube heat transfer performance (compared to the in-line configurations). This may 


indicate bubbles depart TURBO-B tubes differently than smooth tubes. 
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I. INTRODUCTION 


The United States Navy intends to eliminate its dependence 
"ww chlorofluorocarbons (CEGHS) as soon as feasible. 
Currently, the Navy is using CFC-114 (Refrigerant or R-114) in 
most of it's chilled water air conditioning (AC) plants. 
These shipboard AC plants provide the chilled cooling water to 
cool the vital weapons systems and electronic components 
besides providing habitability comfort for the crew. 

The Navy has identified a hydrochlorofluorocarbon, HCFC- 
124 or R-124, as a possible mid-term replacement for R-114 in 
naval shipboard AC systems. With an ozone depletion potential 
(ODP) of only 2$ compared to 70% for R-114, R-124 is a much 
more tolerable fluid to use until an acceptable alternative 
refrigerant can be found that has zero ODP. 

Although both refrigerants operate near the same 
evaporating temperature (2.2 °C), the characteristics of R-124 
pucr similar conditions as R-114 are quite different. For 
example, the refrigeration components (i.e., heat exchangers, 
piping, etc.) must be able to withstand: Cher nirgner 
operating pressures of R-124 (11 psig at 2.2 °C compared to 0 
psig for R-114) and (2) the higher vapor pressure after the 
system is secured (in a main engineering space at 100 °F, R- 
4 would reach 66 psig compared to 31 psig for R-114). Given 


the system components can operate safely at these higher 


pressures of R-124, then the hardware modifications to the 


system should be limited to: 


1. Modifying the compressor to maintain rated capacity and 
use a miscible alkylbenzene oil for lubrication. 


2. Reset/replace the system safety devices to account for 
the higher pressures of R-124. 


3. Replace the system seals with neoprene seals, which 
resist deterioration when in contact with R-122 mm 
maintain leak tight integrity. 


Besides hardware modifications, another main concern is 
how does the heat transfer characteristics of  R-124 


compare to R-114. Table 1, from Bertsch [Ref. 1], provides 


TABLE 1. R-114/R-124 PHYSICAL PROPERTIES AT 2.2 °C 


[Property aa [eaa —— 

[esat (kPa) fe cm q 
Cp (Lig) (J/kg-K) 
k(liq) (W/m-K) 0.0701 0.0744 
p(lig) (g/m+s) 

(N/m) 























g 








the physical properties of R-114 and R-124 at a saturation 
temperature? of 2,2 €; As seen from Table 1, R-124 has 


higher latent heat of vaporization, vapor density and 


thermal conductivity than R-114 at 2.2 °C and a considerably 
lower liquid viscosity. Therefore, it is reasonable for one 
to assume that R-124 would have better heat transfer 
persrormance at 2.2 °C than R-114. Bertsch [Ref. 1] evaluated 
and compared the heat transfer performance of R-114 and R-124 
for a single horizontal smooth tube and various enhanced tubes 
using pure refrigerant and refrigerant/oil mixtures. One of 
his conclusions was: 

For smooth and finned tubes, significant improvements in 
heat transfer (>50%) were obtained in pure R-124 compared 
with pure R-114 at the same saturation temperature. This 
was attributed to the higher saturation pressure of R-124 
activating more nucleation sites and the higher latent 
heat of vaporization, vapor density and thermal 
sensuctrvity of R-124 over R-114. 

Perry [Ref. 2] investigated the effect of oil, pressurization 
and subcooling on the onset of nucleate boiling of R-124 from 
a single horizontal tube. When Perry [Ref. 2] compared his 
results to Bertsch [Ref. 1] for pure R-124, he found a larger 
incipient wall superheat was being obtained, but almost 
identical data trends were observed at high heat fluxes and 
decreasing heat fluxes. He attributed the larger incipience 
to his degassing procedure of the refrigerent system after 
reassembly. 

Now, what if multiple tubes are arranged in a bundle 

configuration? How do the lower tubes influence the heat 
transfer performance of the upper tubes? Lake [Ref. 3] in 


search of an answer to these questions, used a configuration 


of two horizontal tubes (smooth and enhanced tubes) in 
vertical alignment, submerged in a pool of R>114, separa TEE 
a pitch-to-diameter (P/D) ratio. After using various 
combinations of P/D ratios and lower tube heat fluxes, Lake 
[Ref. 3] found that a nucleating lower tube siçnifis s, 
enhanced the heat transfer performance of an upper tube, 
particularly at low heat fluxes. Also, a P/D ratio of 1.8 
gave the best heat transfer performance for smooth tubes at 
high heat fluxes (no P/D effect was found using the enhanced 
tubes). 

To fully understand the nucleate peor pen r g 
characteristics of R-124 and to directly compare the heat 
transfer performance of R-124 to R-114, multiple tube pool 
boiling data for R-124 is needed. Therefore, the following 
thesis objectives were established: 

l. Collect natural convection and pool boiling heat tram s m 
data of pure R-124 using a single smooth and enhaheea 
tube and check the repeatability with the existing R-124 
data of Bertsch [Ref. 1] and Perry [Ref. 2]. 

2. Collect natural convection and pool boiling heat, Amame 
data of R-124 over a range of increasing and decreasing 
heat fluxes using two tubes while varying the spacing 


between the tubes and the heat flux on the lower tube. 


3. Compare the R-124 data to the R-114 data obtained by Lake 
(ne rows 


II. LITERATURE SURVEY 


A. POOL BOILING 

Pool boiling occurs when a heated surface is submerged 
Below the free surface of a liquid pool, at a temperature 
sufficiently above the local saturation temperature to cause 
bubbles to form from the surface. Figure 2.1, from Bar-Cohen 
[Ref. 4], is a pool boiling curve for a highly wetting liquid 
with several regimes displayed. Only the natural convection 
(or single phase convection), boiling incipience and nucleate 
boiling regimes are discussed. 

l. Natural Convection 

The regime (a-b) in Figure 2.1, is where heat transfer 

occurs because of the liquid flow caused by buoyancy forces 
due to density differences in the liquid pool. In this 
regime, the heat flux from the heated surface is proportional 
to some power of wall superheat (i.e., q" « AT” or AT," >”). 
This proportionality remains valid until the wall superheat 
reaches a value where vapor bubbles first form on the heated 
surface. 

2. Boiling Incipience 

The onset of nucleate boiling (ONB) or incipient 


boiling is a random phenomenon (You et al. [Ref. 5]) which 
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Figure 2.1 Pool Boiling Curve 


N-Entiluenced "by factors such as surface condition of the 
heated surface, liquid properties and dissolved gases in the 
mequid pool.  Bar-Cohen and Simon [Ref. 6] and Tong et al. 
[Ref. 7] investigated boiling incipience of highly-wetting 
liquids using refrigerants and dielectrics. They found that 
these liquids effectively flood the micro-cavities on a tube 
surface, causing fewer active embryos (micro-cavitles 
containing trapped vapor). Ultimately, these highly-wetting 
liquids require larger wall superheats to initialize bubble 
growth in the flooded micro-cavities. As a result of the 
larger wall superheats (up to 100 ?C, Bar-Cohen and Simon 
[Ref. 6]), nucleation at (b), in Figure 2.1, occurs in an 
explosive manner. This explosive boiling from the active 
micro-cavities causes che larger inactive cavities to 
activate, rapidly decreasing the required wall superheat to 
a Cin boiling, noted as AT, in Figure 2.1. This drift from 
NEO (cC) only occurs during increasing heat flux and is 
referred to as thermal overshoot or superheat excursion. 
3. Nucleate Boiling 

The regime (c-d) in Figure 2.1 is the fully developed 
nucleate boiling region. At (c), small isolated bubbles are 
leaving the heated surface at low heat fluxes. As the heat 
flux is increased, these small bubbles become more numerous as 
more active sites are created. At high heat fluxes, these 


bubbles begin to coalesce, forming slugs of vapor, which grow 


and depart from the surface. According to Thome [Ref. 8], in 
this regime, heat removal is accomplished by sensible heat in 
the form of superheated liquid or as latent heat in vapor 
bubbles. Also, the principal mechanisms of heat transfer are: 
(1) liquid-phase convection, enhanced by the stripping of the 
thermal boundary layer by the departing vapor bubbles, and (2) 
evaporation of the thin liquid micro-layer trapped between a 


growing vapor bubble and the heated surface. 


B. SINGLE TUBE HEAT TRANSFER EXPERIMENTS WITH R-124 AT NAVAL 
POSTGRADUATE SCHOOL 


Bertsch [Ref. 1] and Perry [Ref. 2] independently studied 
the heat transfer performance of single smooth and enhanced 
tubes submerged in a liquid pool of R-124 at 2.2 °C. Perry 
(Ref, 2] determined that in natural convection the 
experimentally based correlation, equation (1), developed by 
Churchill and Chu [Ref. 9] fairly well predicted the heat 


transfer performance of a smooth tube. 


_ Ra 1/6]° 
Nu « |o. «0 aer E s] | (1) 


where Nu = Nusselt number = hDo/k 


Ra = Rayleigh number = (g:P:(Twall - Tsat) Do?) / (uv-:«) 


Ex nucleate boiling, most correlations available to 
pmaedrot the heat transfer performance of a smooth tube are 
expressed as heat flux as a function of tube wall superheat. 
E het 1) used a best fit) correlation through his 
experimental nucleate boiling data to express the heat flux 


from a smooth single tube in pure R-124 as: 


qu E Um ama (2) 


In the above equation, the wall superheat AT, is expressed in 
degrees Celsius. 

For a single enhanced tube, Bertsch [Ref. 1] concluded 
that nucleate pool boiling enhancement (i.e., improvements in 
heat transfer compared to a single smooth tube) were as much 
as 6.2 times for re-entrant cavity surfaces (TURBO-B and HIGH 
FLUX). Thome [Ref. 10] attributes these enhancements to the 
numerous re-entrant cavities that trap vapor, which create a 
large boiling site density, reducing the wall superheat 


required to sustain nucleation. 


C. SMALL TUBE ARRAY HEAT TRANSFER PERFORMANCE 

The effect of a lower tube on the heat transfer 
performance of an upper tube is very important when designing 
flooded evaporators. It is essential to understand the 


mechanisms that affect the heat transfer performance in small 


tube arrays in natural convection and nucleate boiling for 
both smooth and enhanced surfaces. 
l. Natural Convection 

In natural convection, the interaction effects of a 
lower heated tube on the heat transfer performance of an upper 
tube are the result of two offsetting mechanisms. With a 
lower tube heated, a rising heated buoyant plume impinges on 
and flows around an upper tube, which subjects the upper tube 
to local forced convection effects. This forced convection 
increases the local fluid velocity around the upper tube and 
raises the upper tube heat transfer coefficient: 
Additionally, the heated plume from the lower tube increases 
the local ambient temperature of the fluid impinging on the 
upper tube. This reduces the natural convection buoyancy 
driving force for this tube, decreasing the upper tube heat 
transfer coefficient. Therefore, the presence of a lower 
heated tube has opposing influences on the heat transfer 
performance of an upper tube. Which of these influences is 
dominant depends on the tube spacing, local heat flux and the 
thermophysical properties of the fluid. 

Sparrow and Niethammer [Ref. 11] studied the effect of 
vertical separation between two horizontal cylinders on the 
heat transfer performance of an upper cylinder in natural 


convection using air. It was discovered that heat transfer 
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enhancement! increased as the distance between the tubes 
increased with a maximum enhancement of approximately 1.3 for 
ER P/D< 9: Also, for a P/D < 3, the enhancement decreases 
below 1.0. Al-Alusi and Bushnell [Ref. 12] and Tokura et al. 
[Ref. 13] also working with a small array of horizontal smooth 
cylinders in air, found generally the same results. As the 
distance between the cylinders increased, the heat transfer 
enhancement increased, reaching a maximum, then returning to 
a lower value of enhancement. 

Sparrow and Boessneck [Ref. 14] investigated the 
effect of transverse misalignment of a lower cylinder on the 
heat transfer performance of an upper cylinder, for 2 s P/D s 
ONE Sing air. For the perfectly aligned upper and lower 
cylinders, the heat transfer enhancement of the upper cylinder 
increased as the P/D ratio increased. Again, the largest 
enhancement was found with a P/D of 9 and the P/D of 2 
produced the lowest  enhancment. As  misalignment was 
increased, the enhancement began to merge to a common value of 
about 1.03 at sufficiently large offsets (about 2.5 cylinder 
diameters). Sparrow and Boessneck concluded that this 
behavior is caused by the rising plume from the lower cylinder 


presenting a pre-heated quasi-forced-convection flow to the 


lEnhancement is defined for a small tube array as the 
upper cylinder heat transfer coefficient when both cylinders 
are heated divided by the upper cylinder heat transfer 
coefficient when the lower tube is unheated. 
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upper cylinder. At smaller P/D ratios, the increased local 
ambient temperatures of the fluid impinging on the upper 
cylinder dominates and decreases the upper cylinder heat 
transfer performance. At larger P/D ratios, the Momeed 
convection effects are more dominant and increase the heat 
transfer from the upper cylinder. 

Inagaki and Komori [Ref. 15] studied natural 
convection heat transfer from two vertically aligned 
horizontal cylinders in water. As discovered with air, a 
decrease in enhancement which falls below 1.0, is found for 
P/D < 1.0. A maximum enhancement of 1.25 occurs between 2.5 
< P/D s 3.0 and at larger P/D ratios, the enhancement 
decreases to 1.0. Inagaki and Komori attribute this behavior 
to the same two opposing influences discussed previously. 

2. Nucleate Boiling 

For small tube arrays in nucleate boiling, the most 
Significant mechanism of enhancement within the array, is the 
well documented bundle effect”. Cornwell and associates 
(Refs. 16-19] studied fluid flow in tube bundles and found 
Significiant enhancements in heat transfer in the upper 
portions of the bundles, concluding that forced convection and 


bubbly flow from the lower tubes, creating sliding bubbles 


“Bundle effect is the enhancement of an upper tube due to 
the activation (i.e., the nucleation of bubbles) of lower 
tubes within the bundle. 
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around the upper tube, cause significant enhancement at low 
heat flux conditions. 

Fujita et al. [Ref. 20] conducted tests with an array 
E three horizontal smooth tubes in R-113 at 0.1 MPa. For 
heat fluxes « 20 kW/m’, the bundle effect enhanced the upper 
tube heat transfer performance. Fujita et al. [Ref. 20] also 
used an array of two smooth tubes, where the lower tube heat 
flux was held constant and the upper tube heat flux was 
decreased uniformly. The upper tube enhancement increased as 
the lower tube heat flux was raised. Lake [Ref. 3] 
investigated the effect of P/D ratio and lower tube heat flux 
on the heat transfer performance of an upper tube, using an 
array of two smooth tubes in R-114 at 1 atmosphere. For these 
tests, the bottom tube heat flux was held constant and the 
upper tube heat flux was raised and lowered steadily. Lake 
found that for an upper tube heat flux of « 20 kW/m/, the 
enhancement of the upper tube was due to a bundle effect and 
the amount of enhancement increased with increased lower tube 
heat flux, which is consistent with Fujita et al. [Ref. 20]. 
Also, Lake found that for upper tube heat fluxes of > 20 
W/m, all bundle effects faded and that the P/D ratio of 1.8 
provided the best upper tube heat transfer performance. 

For a small array of enhanced tubes, several authors 
report findings similar to those found for the smooth tube 
arrays. Stephan and Mitrovic [Ref. 21], using a small array 


of three horizontal T-shaped finned tubes in R-114 at 0.2 MPa 
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and R-12 at 0.15 MPa, found bundle effects present at heat 
fluxes < 20 kW/m’, with these effects disappearing as the heat 
flux was raised. Hahne and associates [Refs. 22-24], \eame 
two different fin arrays in R-ll at 0.1 MPa, also me mma 
significant bundle effects for heat fluxes of « 20 kW/m, 
which subside at higher heat fluxes. Hahne et al. [Ref. 24] 
also investigated the influence of P/D ratio on upper tube 
heat transfer. For heat fluxes « 20 kW/m^, the maximum P/D 
ratio tested (P/D of 3.0) produced the best heat transfer 
performance for both finned arrays. Above 20 kW/m?, neither 
P/D ratio nor tube fin pitch affected the upper tube heat 
transfer. Lake [Ref. 3], in conjunction with his testing of 
smooth tube performance above, found no effect of P/D ratio on 
HIGH FLUX tube (porous coated tube) performance at any heat 
flux and that there was no effect from a lower heated tube on 
the heat transfer performance of a nucleating upper HIGH FLUX 
tube. Fujita et al. [Ref. 20], like Lake, also tested porous 
coated tubes using the same testing conditions as for the 
smooth tubes and he concluded that the lower tube heat flux 
setting had no influence on the heat transfer coefficient of 
the upper tube, which agrees with the conclusions of Lake 


[Ref. 3]. 
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III. EXPERIMENTAL APPARATUS 


A. SYSTEM DESCRIPTION 

A schematic diagram of the experimental apparatus is shown 
in Figure 3.1. Details of the apparatus design, construction 
and calibration are discussed by Karasabun [Ref. 25], Reilley 
[Ref. 26] and Sugiyama [Ref. 27]. The apparatus is very 
Sar Co the one used] by Lake [Ref. 3] to conduct pool 
boiling of pure R-114 with a two tube configuration. The only 
major change being the replacement of the Pyrex-glass Tees 
used as the condenser and evaporator. Bertsch [Ref. 1] 
replaced the Pyrex-glass Tees with fiberglass reinforced 
Pyrex-glass Tees enabling the system to withstand the 
increased operating pressure of R-124. 

As shown in Figure 3.1, the present system contains the 
following components: (1) a reinforced Pyrex-glass Tee, used 
as the evaporator to boil the refrigerant; (2) a reinforced 


Pyrex-glass Tee, used as the condenser to condense the 


refrigerant vapor; (Sia liquid serrigerant | storage 
reservoir; (4) an oil storage and mixing subsystem (not 
used); (5 a vacuum pump; (6) an ethylene-glycol/water 
mixture cooling subsystem; (7) a condenser pump; 


(8) aluminum framework with plexiglass siding which houses 
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Figure 3.1 Schematic of Experimental Apparatus 
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Sempoments (1) through (4); and (9) a data acquisition and 
IUSErumentation system (not shown in Figure 3.1). 

The system operates in a closed loop with vapor being 
generated by the heated tubes submerged in the R-124 liquid 
pool in the evaporator. The vapor travels to the condenser, 
via a two inch diameter aluminum L-shaped pipe, where it's 
condensed on a copper coil cooled by the chilled ethylene- 
glycol/water mixture. The condensed liquid refrigerant is 


gravity fed back to the evaporator via the return line. 


B.  BOILING TEST SECTION 
l.  Evaporator 

The evaporator is a Corning Pyrex-glass Tee section 
(three inch nominal inner diameter) coated with a continuous 
winding of fiberglass impregnated filaments bonded with a 
modified epoxy resin. It was installed horizontally with the 
Side-arm of the Tee vertical and hydrostatically safety tested 
Boo psig by Bertsch [Ref. 1]. Aluminum endplates, 210 mm 
diameter and 12.7 mm thick, house the thermocouple wells, an 
oil entry fitting (not used) and the Teflon inserts, which 
hold the test tubes in position. The endplates are mated to 
the evaporator using cast-iron flanges. At the very bottom of 
the evaporator sits the refrigerant liquid return line from 
the condenser. A copper deflector plate was installed just 


above the return line to prevent any vapor bubbles that might 
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rise off the return line from impinging on the test tubes. The 
liquid refrigerant level in the evaporator is approximately 10 
mm below the top of the evaporator and approximately 10 mm 
above the top of the upper tube. Figure 3.2 is a sketch of 
the evaporator and a cast-iron flange. 
2.  Evaporator Tubes 

Figure 3.3 is a schematic of a boiling kep Two 
different tube types were used in this study: a smooth tube 
and a TURBO-B enhanced surface tube. Table 2 provides the 


dimensions of each tube type while Figure 3.4 shows an 


TABLE 2.  DIMENSIONS OF TEST TUBES 





mm -— EP —À a  J——————————— ——— 


E ICH uL 
Pruaso-8[11.6]13.8)11.8 | 15.9] 203.2] 76.2| 299 | 


| TURBO-B | 11.6 | 13.8 |11.8]15.9]| 203.2| 76.2] 398 





expanded view of the enhanced surface of the TURBO-B tube. 

Each tube is heated by a 1 kW, 240 volt stainless-steel 
Cartridge heater. The heater measures 6.35 mm outside 
diameter, 203.2 mm overall length with a heated length of 190 
mm and is centered inside the test tube. The tubes are 


instrumented with eight thermocouples at various locations to 


18 





Corning Pyrex Glass Evaporator (D x N = 402x51 mm 
L = 178 mm, Li= 127 mm) 


, 





Cast Iron Flange (d, = 190 mm, d, = 210 mm) 


Figure 3.2 Sketch of Pyrex Glass Evaporator 
and Cast-Iron Flange 
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monitor the tube wall temperatures. Bertsch [Ref. 1] provides 
a detailed assembly procedure for the manufacturing and 
instrumentation of a test tube. 

The test tubes are secured in the evaporator by Teflon 
inserts and sealed with o-rings (Figure 3.5). Strong-backs, 
shown in Figure 3.6, are used to anchor the Teflon inserts in 
the endplates. 

3. Tube Configurations 

Figure 3.7 shows the six tube configurations that were 
used. Configuration CO is the single tube configuration used 
by Bertsch [Ref. 1] and Perry [Ref. 2]. In configuration! 
through C4, the top tube is fixed relative to the evaporator 
while the lower tube position is varied. Configuration C5 
consists of rotating the Teflon insert used for C4 through an 
angle of 30°. Note that when rotation is completed, the upper 
tube position has moved off the evaporator centerline and the 


lower tube position remains unchanged. 


C.  CONDENSER SECTION 

The condenser is also a fiberglass-epoxy coated Corning 
Pyrex-glass Tee section identical in size, shape and substance 
to the evaporator show iqure p n It is mounted 
vertically (see Figure 3.1) with aluminum endplates bolted to 


cast-iron flanges, similar to the evaporator arrangement. The 
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Figure 3.5 Sketch of an Endplate and Teflon Insert 


23 






TTT 


LIN 


` B 
sË 





suor3eanbrguoo eqnj L°¢ sanbry 


(qa55nid) g38nL ON - O “HAND - 6 


(,0€ daLbvbou ‘z2°t d/d) (c't a/d) 
Go po 
d 
(8't a/d) (o'z a/a) 
23 TO 


ONO 


LE 


(s`t q/d) 
£O 





02 


25 


top endplate is fitted with connections for refrigerant 
charging/vacuum operations and the ethylene-glycol/water 
coolant inlet tube. The bottom endplate is fitted with 
connections for the gravity feed return line to the evaporator 
and the ethylene-glycol/water coolant outlet tube. The 
refrigerant vapor rises from the evaporator, passing through 
the L-shaped aluminum pipe and enters the condenser where it 
is condensed on a 3/8 inch outside diameter copper tube which 


is helically coiled inside the condenser. 


D. REFRIGERANT STORAGE RESERVOIR 

The refrigerant storage reservoir is an aluminum cylinder, 
230 mm in diameter and 254 mm in height. It is located within 
the apparatus enclosure above the evaporator and below the 
condenser (see Figure 3.1). The reservoir enables the 
operator to store the liquid R-124 instead of releasing it 
into the atmosphere, prior to changing out the test tubes. 
The installed sight glass allows monitoring of the liquid 
level during replenishing from an outside source (i.e., 145 
lb. storage cylinder) or refilling the evaporator after atest 
tube changeout. The reservoir is isolated from or aligned to 


the system by means of valves V6, V7, V12 and V13. 
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E. OIL SUBSECTION 

This study was devoted to pure R-124. In order to ensure 
that no oil contaminated the system, the oil line from the oil 
cylinder to the evaporator was removed (see Figure 3.1) and 
the access on the evaporator endplate was capped. Perry 
[Ref. 2] describes the oil subsection components and the oil 


addition procedure for a refrigerant/oil mixture study. 


F. COOLING SUBSYSTEM 
1. Ethylene-Glycol/Water Coolant Sump 
At the base of the apparatus sits a 0.15 cubic meter 
plexiglass tank containing approximately 30 gallons of 
ethylene-glycol/water mixture producing a freezing point of 
-25°C. The sump rests on a wood platform and is wrapped with 
22 mm thick foam rubber sheet for insulation. Supply and 
discharge connections for the condenser and sump pump as well 
as a thermocouple to monitor the sump temperature are 
installed. 
2. Refrigeration Plants 
A 1/2 ton R-502 and 1/4 ton R-12 refrigeration plants 
are used simultaneously to maintain the ethylene-glycol/water 


coolant temperature at about -15? C. Both plants are located 


close to the sump and contain the following components: (1) 
hermetically sealed compressor assembly; (2) air-cooled 
condenser; (3) receiver; (4) filter-dryer; (5) pressure 


2 


regulator; (6) temperature control switch; and EE) 
thermostatic expansion valve. The R-502 plant (Figure 3.8) 
uses a counter flow heat exchanger as the system evaporator 
while the R-12 plant uses a 9.5 mm diameter coiled copper 
tubing submerged in the sump tank as the system evaporator 
(Figure nl M 

3. Pumps and Control Valve (VC) 

A positive displacement sump pump, rated at 8 gpm, is 
connected to the discharge connection of the coolant sump with 
a one inch diameter PVC pipe. It circulates coolant from the 
sump through the R-502 evaporator and back to the sump. A 
similar pump, called the condenser pump, circulates coolant 
from the sump to the condenser, via the control valve VC,and 
back to the sump (see Figure 3.1). Valve VC is used by the 
operator to control the refrigerant pool saturation 
temperature. A manually operated bypass valve (V9, Figure 
3.l) is installed upstream of valve VC to prevent overloading 


the condenser pump when valve VC is closed. 


G. APPARATUS ENCLOSURE 

A welded aluminum frame measuring 214 cm high, 61 cm wide 
and 5l cm deep is used to contain most of the system 
components. The upper half of the frame contains the 
evaporator, condenser, oil subsystem, and refrigerant storage 


reservoir (see Figure 3.1). The lower half of the frame 
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contains the coolant sump tank. For operator safety and to 
reduce heat leakage into the evaporator and condenser, all 
sides of the upper half of the frame are covered with 13 mm 
thick plexiglass, hinged on both sides to allow access to the 
system components. Holes in the front plexiglass cover allow 
operation of valves V1 through V8 while maintaining the 


integrity of the enclosed upper frame. 


H. INSTRUMENTATION 
l. Power Measurement System 
Figure 3.9 is a schematic of the power measurement 
system. Each boiling tube heater is powered by a 240 volt AC 
source controlled by a variac. The variac output ranges from 
0-220 volt and 0-5 amp, adjustable by the operator to obtain 
the desired heat flux at the tube surface. The current and 
voltage sensors convert the AC input signal into DC voltages 
which are used as input by the Data Acquisition Unit. 
2. Temperature Measurement 
Temperatures are measured using 30 gauge (0.25 mm 
diameter) copper-constantan (type T) thermocouples. For each 
test tube, eight thermocouples are spaced over the heater 
sleeve, as shown in Figure 3.10, to measure the wall 
temperatures. Three thermocouples are inserted into housings 


that penetrate the evaporator endplate and are submerged in 
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the refrigerant liquid pool at different locations (see Figure 
ESI. Each housing is constructed of a low conductivity 
Epcundess steel bedy and a Nigh conductivity copper tip 
(Figure 3.11). This minimizes the errors in pool temperature 
Measurements caused by heat conduction through the housing. 
FR single thermocouple is located in the coolant sump tank to 
measure the sump temperature. 

Thermocouples were read by the a Hewlett-Packard (HP) 
3497A Data Acquisition Unit which is controlled by an HP 326 


Computer System. 


1. DATA ACQUISITION AND REDUCTION PROGRAM 
Data acquisition and reduction program DRP72, used by Lake 
[Ref. 3], was modified to incorporate the thermophysical 
properties of saturated R-124. Bertsch [Ref. 1] discusses the 
generation of the equations used in the program calculations. 
The modified DRP72 program was renamed DRPJY. Appendix B 
contains the program DRPJY. The program DRPJY directs the HP 
362 Computer System and the HP 3497A Data Acquisition Unit to 
read and process all thermocouple, current and R.M.S. voltage 
EEG ngs. Using this information, DRPJY performs the 
following calculations (from Bertsch [Ref. 1]): 
imeconverts voltage readings into tube temperatures and 
power. 


2. Computes heat transfer rate from cartridge heaters. 
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3. Computes average test tube wall temperatures and wall 
superheat (Twall - Tsat). 


4. Computes physical properties of R-124 using property 
correlations from  Reteror. r [Ref- 28] at the film 
temperature (Twall + Tsat)/2. 


E Computes heat transfer coefficient of R-124 from the 
unheated ends and heated length of the test tubes. 


6. Computes heat losses from the unheated ends of the test 
IUOS. 


7. Computes heat flux from heated length of test tubes. 


The processed output data was printed by an HP Inkjet 
printer and stored on the computer hard drive. As a 
precaution, each output data file was copied to a 3.5 inch 
floppy disk in the event the computer hard drive becomes 
unreadable. Table 3 lists the channel assignments that are 


sensor inputs to the Data Acquisition Unit. 
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TABLE 3.  HP-3497A CHANNEL ASSIGNMENTS 


DOLO Upper test tube wall temperature (T1-T8) 
SE Lower test tube wall temperature (T9-T16) 
Liguidip poo Ethe rmocoup k 


Channel Assignments 
08 
) 
Liquid pool thermocouple (Tld2) 


Liquid pool thermocouple (iv 


| 15 | Liquid pool thermocouple (Tv) — 
DUNE COO EMPERO 


21 R.M.S. voltage of tube heaters | 


Upper tube current sensor 
22 


Lower tube current sensor 





36 


A. 


IV. EXPERIMENTAL PROCEDURE 


EVAPORATOR DISASSEMBLY PROCEDURE 


The following procedure is used to transfer the liquid 


refrigerant from the evaporator to the reservoir to facilitate 


test tube changeouts: 


I. 


Operate both R-502 and R-12 refrigeration systems to 
cool the sump tank temperature to approximately -15 °C. 


Start the condenser pump and using control valve VC, 
adjust the coolant supply pressure to 1 psig. 


Close valve Vô to isolate the evaporator from the 
condensing refrigerant vapor (see Figure 3.1). 


Align the refrigerant storage reservoir to the system by 
opening valves V6, V7, V12 and V13. 


Position and energize the high volume air dryer (set for 
high heat) so the heated air stream is blowing on the 
center of the evaporator. Adjust valve VC to ensure the 
refrigerant vapor is condensing. 

After 1-2 hours, the evaporator should be dry with the 
liquid refrigerant in the reservoir and a small amount 
in the return line. 


Close valves V6, V7, V12 and V13, isolating the 
reservoir from the system. 


Secure and remove the air dryer. 

Stop the condenser pump. 

Open valves V8 and V10 to reduce to system pressure to 
atmospheric. Close valve V8 and V10 when the system 


pressure reads 0 psig on the system pressure gage. 


Disconnect the test tube heater and thermocouple leads. 
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12. Remove strongbacks. 


13. Remove Teflon inserts and test tubes. 


B. EVAPORATOR ASSEMBLY PROCEDURE 
1. Evaporator Assembly 

Before installing smooth test tubes, rinse with 
acetone then wipe clean with lint-free wipes or if installing 
TURBO-B test tubes, brush the enhanced surfaces with a soft 
bristle brush (i.e., toothbrush) using ethyl alcoholiithen 
repeat the brushing using acetone and let the tubes aims 
Check the Teflon insert o-rings for damage and replace as 
necessary. Before installing the test tubes, wipe the 
evaporator interior clean with lint-free wipes. Install the 
test tubes in the evaporator using the Teflon inserts to 
secure the tubes in place. Snug the strongbacks up to the 
Teflon inserts and connect the thermocouple leads. Prior to 
connecting the heater leads, check the resistance between the 
heater leads and the outer tube wall using a multimeter. One 
megohm resistance or greater is satisfactory. 

2. Leak Test 

After assembling the evaporator and before charging 
the system with refrigerant, perform a vacuum test to verify 
the integrity of the system. Evacuate the system by starting 
the vacuum pump and opening valves V11 and V8 (see Figure 


3.1). After evacuating the system to approximately 29 
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IO, close valves Vll and Vê, “then stop the vacuum pump. 
Beene drop 1n vacuum is detected (i.e., mo leak) within 30 
memaces, Charge the system with R-124. 

If a leak is detected, open valves V8 and V11 and 
Pe@essurize the system with air to 20 psig. Apply a soap/water 
solution to the system sealing surfaces to identify the leak 
location (bubbles will be seen from the leak source). After 
identifying and repairing the leak, retest the system for 
leaks. 

3. Charging System with Refrigerant 

Once the system integrity is assured, fill the 
evaporator with R-124, using one of the following procedures, 
to a level just below the top of the evaporator (NOTE: some 
refrigerant will be lost during the degas procedure, 
Buenetore, it is important to fill the evaporator above the 
level needed to conduct data runs): 

a. From Storage Reservoir 

Open rer EES 
between the storage reservoir and the system. Open valves 
ieee vo and V5 to gravity fill the evaporator with liquid 
refrigerant. When the evaporator is filled to the proper 
level, close valves V6, V13, V7 and V12 to isolate the storage 


reservoir from the system. 


59 


b. From Refrigerant Cylinder 

Connect the cylinder to valve Vl0N usrnomm 
charging hose (Figure 3.1). Open valves V5, V8 and NUMEN 
then open the cylinder liquid valve to release the refrigerant 
from the cylinder into the system. Close valves V8, V10 and 
the cylinder liquid valve when the evaporator is filled to the 
proper level. Disconnect the charging hose from valve V10 and 

stow the cylinder. 

4. Degassing Procedure 

This degas procedure is similar to Perry's [Ref. 3] 
degassing procedure B. Degas the system by boiling the 
refrigerant, using a upper test tube heat flux of 85 KW/m^ and 
a lower test tube heat flux of 20 KW/m, for 30s»minutes: 
Ensure the condenser pump is running and control valve VC is 
Open to prevent overpressurizing the system. During the 
boiling process, air and non-condensible gases rise in the 
system and collect in the condenser. After the 30 minute 
boiling time has elapsed, secure the heat flux to the test 
tubes and stop the condenser pump. Let the system set for 30 
minutes allowing the condensed refrigerant to drain to the 
evaporator. After the 30 minute set time has elapsed, start 
the vacuum pump and open valve Vll, then open valve V8 for 20 
seconds to remove the accumulated air and non-condensible 
gases. After closing valve V8, stop the vacuum pump and close 


valve V11. 
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5. Data Acquisition Channel Check 
After all evaporator thermocouples (test tubes and 
pare connected, load and run program SETUPJY to verify 
the output of each thermocouple is correct. Any output that 
is inconsistent or suspect is either fixed or noted and later 


excluded from the calculations in program DRPJY. 


C. STANDARD OPERATING PROCEDURE 
The following procedure is used to conduct each data 


zo lection run: 


1. Operate the R-502 and R-12 refrigeration systems to cool 
the sump to approximately -15 °C. 


2. Energize the HP-362 Computer System, Data Acquisition 
Unit and the variac control panel. 


3. Load and run program SETUPJY. Verify the thermocouple 
GUEDUES are consistent. 


4. Verify control valve VC is closed, then start the 
condenser pump. 


5. Slowly open valve VC until the coolant pressure gage 
reads 1/8 psig. This pressure setting will ensure that 
the refrigerant pool will be slowly cooled. 


6. Monitor the thermocouple outputs with SETUPJY until 
confident the thermocouples are operating properly. 
Note which thermocouples are inconsistent and must be 
excluded from program DRPJY. 


7. When the refrigerant pool temperature is approaching 5 
°C, load and run program DRPJY. Select 'TAKE DATA' 
ept10n. 


8. At the prompt, enter the filename that the data is to be 


stored under (page 38 explains the file naming scheme 
used). 
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LS 


duis 


P. 


le. 
14. 


O 


O 


NE 


ach 


OR 


D. 


At the prompt, enter the number of defective upper tube 
thermocouples that were identified by SETUPJY. 


At the prompt, enter the defective upper tube 
thermocouple locations that were identified by SETUPJY. 


At the prompt, enter the number of defective lower tube 
thermocouples that were identified by SETUPJY. 


At the prompt, enter the defective lower tube 
thermocouple locations that were identified by SERPENS 


Select the type of test tubes installed. 


From the main menu, select 'SET Tsat'. Monitor the pool 
thermocouples and, if required, adjusted valve VC to 
maintain the desired pool saturation temperature of 2.22 
E AER 


From the main menu, select 'SET UPPER TUBE HEAT FLUX', 
then adjust the main variac to the desired heat flux. 
Initial heat flux setting is 500-600 KW/m’. 


Press the keyboard button fl to return to the main menu 
and select 'SET LOWER TUBE HEAT FLUX'. Adjust the aux 
variac to the desired heat flux. 


From the main menu, select 'SET Tsat'. Adjust valve VC 
to maintain the desired pool saturation temperature 
while monitoring the pool thermocouples. 


Once the pool saturation temperature is stable, return 
to the main menu and select 'TAKE DATA'. 


For each heat flux setting steps 15-18 are repeated. 


DATA FILES 


Data were filed using the following convention: 


Example (D0904ASC1): 


iig 
Zs 


First letter D indicates the file is a data file. 


The following four characters indicate the date; in this 
example September 04, 1993. 
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3. The next letter indicates the order in which the data run 
was completed on that day; no letter - lst run, A - 2nd 
Wn. B jr un, etc. 


4. The next one or two letters indicate the tubes used; 
S - smooth, TB - TURBO-B. 


5. The last letter and number combination represent the tube 
configuration used. 


Table 4 is a summarized list of the data runs completed. 
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TABLE 4. LIST OF DATA RUNS 


à 





Data File 


D0728ASCO 
D0804ASCO Repeated Ley Po u Ce 


DO816SCO0 Repeatabi E singer 

ID US Smooth Repeatability; P/D = 2.0; no 
power on lower tube 

D0714SCl Smooth Repeatability; P/D = 2.0; no 
power on lower tube 

DOZ S CN Smooth Data; repeatability; P/D = 2.0; 
no power on lower tube 

D0901SC1 Smooth Data; P/D = 2.0; 1 kW/m’ on 
lower tube 

D0902SC1 Smooth Data; P/DE-2270/959 Wm 
lower tube 

D0902BSC1 | Smooth Data; P/D = 2.0; 201 kW/m fox 
lower tube 

D0903SC1 Smooth Dara D 0 m4 onbo 
tubes inc/dec together 

DO904SC1 Smooth Data; P/D'=127 0:05 EkW/m Nom 
upper tube 

DO904ASC1 | Smooth Data; P/D = 2.0; 10 kW/m’ on 
upper tube 

D0905SC1 Smooth Data; P/D = 2.0; 20 kW/m Rom 
upper tube 

D0906SC2 Smooth Data; P/D = 1.8; no power on 
lower tube 

D0906ASC2 | Smooth Data; P/D = 1.8; 1 kW/m’? on 
lower tube 

poSsensc» Smooth Data; P/D = 1.8; 5 kW/m* on 
lower tube — 


D0907ASC2 Smooth Data; P/D = 108; 10. kwW/m lem 


Tube Type | Purpose 





Repeatability; single tube 





lower tube 
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TABLE 4. 


| Data File 
D0908SC2 


D0908ASC2 


D0908ASC2 


BO909SC2 


D0909ASC2 


[091 05SC3 


D0910ASC3 


RDSEISCJ 


DOSH SROS 


EODSISASC3 


DO920SC3 


D0920ASC3 


BOSZOBSC3 


BOSZ20E5C3 


DO921SC4 


D0921ASCA4 


D0921BSC4 


LIST OF DATA RUNS (continued) 


Tube Te Purpose 


Smooth Datan P/D = Ja Ei. Jet At" com 
lower tube 

Smooth DATA eege consboth 
tubes inc/dec together 

Smooth Date wap D Das IE TH a 
upper tube 

Smooth Dataa P D- = Ju; 0 d imeson 
"Open tube 

Smooth Dataa AD — dados 20 0 KWAm on 
lower tube 

Smooth Data; P/D = 1.5; no power on 
lower tube 

Smooth patas p/D = JS al um Pon 
lower tube 

Smooth Parar Ped = SS PKR mon 
lower tube 

Smooth Bata e/D = 1252 S08, maon 
lower tube 

Smooth Data:  P/D = UNS Wo DEI mon 
lower tube 

Smooth Data; P/D = 1.8; q" on both 
tubes inc/dec together 

Smooth Data Me/D= SS kW /maomn 
upper une 

Smooth Bat aqp /D d NONI DIET AOT 
upper tube 

Smooth A = les LON imo 
upper tube 

Smooth Data; P/D = 1.2; no power on 
lower tube 

Smooth Datan B/D = Esc kW/m on 
lower tube 

Smooth Data A a on 
lower tube 
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TABLE 4. LIST OF DATA RUNS (continued) 


——snA— ` ———— —  ————————— 


| Data File Tube Type 


D0921CSC4 Smooth Data; "PAD = 1.2 010 kw/m Fom 
lower tube 


D0922SC4 Smooth Data; P/D = 1225870 EW/m ES 
lower tube 

D0922ASC4 Smooth Data; P/D - 1:27 q 7 on potim 
tubes inc/dec together 

D0922BSC4 Smooth Data: P/D =EN Wn KN, m EO 
upper tube 

D0922CSC4 Smooth Data; P/D —-201 72798 0D KW məm 
upper tube 

D0922LSC4 | Smooth Data; "P/Dp = 172 0PkW/m Mom 
upper tube | 

D0930SC5 Smooth Data; "P/D = 1792" rotated=30 TEE 
power on lower tube 

D0930ASC5 Smooth Data;"P/D -9M2Wrotatedw 0 ia 
kW/m^ on lower tube 

D0930BSC5 Smooth Data; "P/D = 1*?2NFrotatecWNon ES 
kW/m^ on lower tube 

D0930CSC5 Smooth Data; P/D = 12 rotated 0 TENE 
kW/m^ on lower tube 

D0929CSC5 Smooth Data;$P/D = 1 2Protatedes UT 
kW/m^ on lower tube 

D0929DSC5 | Smooth Data; P/D = 1.2 rotated 30°; q" 
on both tubes inc/dec together 

D0929BSC5 Smooth Data; P/D = 1.29" rotated 30 
kW/m’ on upper tube 

D0929ASC5 Smooth pata PP/D = IT Erota ede WE 
kW/m^ on upper tube 

D0929SC5 Smooth Data, aP D A EUME Oot a CUMU 

| kW/m*^ on upper tube 


D11012 p5g0 TURBO-B Repeatability, single tube 
DIOTZATBEO TURBO. Repea tab? bie, Sing estas 
DIQIZBTBCO | TUREO-B Repeatability; single tube 
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TABLE 4. LIST OF DATA RUNS (continued) 





Bn, Bis File L | 
| D1019TBCl 
pr 
D1018CTBCl 
DHUOISDTBCI 
D1019TBCl 
D1019ATBCl 
D1019ETBCl 
D1019CTBC1 
BRUM SIDTBCA 
PRO LOBTBCI 
BROZUTECZ 
D1020ATBC2 


IMM)ZÜBIBC2 | TURBO-B 


== HH 


IMADZOCTBGZ | TURBO-B 


D1020DTBC2 | TURBO-B 


Purpose 


Repeatability; P/D = 2.0; no 
power on lower tube 


Repeatability; P/D = 2.0; no 
power on lower tube 


Data; repeatability; P/D = 2.0; 
no power on lower tube 


Darafrp/D = 02. 05» 1 kumon 
lower tube 


DE P/D == 2005850 KW /mmon 
lower tube 


Pata A2 Oros Amason 
lower tube 


Baran =D don mon 
lower tube 


Data /D = 2400. q" aonsboth 
tubes decreased together 


a A ADA AO 
upper tube 


Data aa Di 92-0] oa mom 
upper tube 


Data? P/D = 2.0:920X9kwW/m^€*on 
upper tube 


Data; P/D = 1.8; no power on 
lower tube 


Data AS co PINK mon 
lower tube 


Data P 45 — ole > a, on 
lower tube 


Data P/D mI B Orki m on 
lower tube 


Parade <il Omk mcm 
lower tube 


TABLE 4. LIST OF DATA RUNS (continued) 


Data File Purpose 

D1020ETBC2 | TURBO-B Data; P/D = 1.8; "omes eM 
tubes decreased together 

DOF vies CZ TURBO-B DS E 10.0; no KM MO 
upper tube 

D1021ATBC2 | TURBO-B Data; P/D = 1.8; 10 kw/m” on 
upper tube 

D1021ATBC2 | TURBO-B Data; P/D = 1.8; 20 kW/m” on 
upper tube 

DIOZTCTBC3) tURee-B Data; P/D = 1.5; no power on 
lower tube 


D1021BTBC3 | TURBO-B Data; P/D - ES; kW mo 
lower tube 


D1021ETBC3 | TURBO-B Data 8P/D = MIS m Mon 
lower tube 


D1021FTBC3 | TURBO-B Data; *PP/D = W25-TOMkW mom 
lower Rune 


D1022TBC3 | TURBO-B Data; P/D = M. 57m2 0EkW m EON 
lower tube 

D1022ATBC3 | TURBO-B Data; P/D = WT 75; 0 MBOnN DO 
tubes decreased together 


D1022BTBC3 | TURBO-B Data; P/D = 15>) 5ekW/m Nor 
upper tube 

D1022CTBC3 | TURBO-B Data; P/D 2-715529 0NEW mmc 
upper tube 


D1022DTBC3 | TURBO-B Data” P/D = q€*5529» 0MkW/ mom 
upper tube 
D1022ETBC4 | TURBO-B Data? P/D = 1.2; mo power on 
lower tube 
D1022FTBC4 | TURBO-B Dat P/D =- M2 Wim Mon 
I lower tube 
D1025TBC4 | TURBO-B {pata age = 1 O° omen 
lower tube 
D1025ATBC4 | TURBO-B Data; P/D.- 1 2:00 KW miem 
lower tube 


TABLE 4. LIST OF DATA RUNS (continued) 


Data File | Tube Type | Purpose | 


D1025BTBC4 | TURBO-B Datan PD -= 1.2 20 kW/m on | 
lower tube 


pAID25CTBC4 | TURBO-B Data PD =". q” on both 
tubes decreased together 

D1025DTBC4 | TURBO-B Data; P/D = 1.2; 5 kW/m on 
upper tube 

D1025ETBC4 | TURBO-B Data; P/D = 1.2; 10 kW/m* on 
upper tube 


D1025FTBC4 | TURBO-B Data: PDS PO Oui mom 
upper tube 


D1026TBC5 TURBO-B Dara; P/D- 12 Trotrtated 1309 mno 
power on lower tube 


D1026ATBC5 | TURBO-B Data pp otra emo dl 
kW/m^ on lower tube 


D1026BTBC5 | TURBO-B Dita MD A MOL tedio 5 
kW/m^ on lower tube 

D1026CTBC5 | TURBO-B Data p Oe J Jy aro pe S p. hp 
kW/m^ on lower tube 

D1026DTBC5 | TURBO-B Data. e5/ De— de carotases ROCA 20 
kW/m^ on lower tube 

D1026ETBC5 | TURBO-B Data; P/D = 1.2 rotated 30°; q" 
on both tubes decreased 
together 

D1027TBC5 TURBO-B Dat? IDE EE 
kW/m? on upper tube 

D1027ATBC5 | TURBO-B Dara DO rorared Bo am 
kW/m? on upper tube 


D1027BTBC5 | TURBO-B EPA cca NO eri teta o 
kW/m? on upper tube 











m 
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V. DISCUSSION OF RESULTS 


A.  REPRODUCIBILITY/REPEATABILITY 

To verify that the equipment and operating procedure being 
used were consistent with those used by Bertsch [Ref. 1] and 
Perry [Ref. 2], six reproducibility/repeatability runs were 
performed. Program DRPGP (from Perry's single tube work 
[Ref. 2]) was used for the runs conducted with a single test 
tube and DRPJY was used for all runs conducted with two test 
tubes. Figures 5.1 and 5.2 are plots of increasing and 
decreasing heat flux using a smooth tube in configuration CO 
(which is identical to the configuration used by Perry 
[Ref. 2]). In Figure 5.1, the three runs completed are in 
very good agreement with each other. At lower heat fluxes 
( < 4 kW/m’) some scatter is present which disappears at mid- 
range heat fluxes (4-8 kW/m’). Perry [Ref. 2] found the same 
phenomenon and attributed the scatter to convection regime 
thermal gradients where mixing is inefficient within the 
liquid pool, after discovering large local variations in the 
tube wall and pool thermocouples of up to 5$. In this Su: 
local variations in the tube wall and pool thermocouples were 
found to be up to 7% at lower heat fluxes and decreasing as 


heat flux increased, which supports Perry's theory. In 
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Figure 5.1 Increasing Heat Flux Repeatability Runs using 
Smooth Tube Configuration C0 


Su 


natural convection, the three data runs started to the left of 
the Churchill and Chu correlation (Ref. 9] “shown as ENE 
line) at low heat fluxes and converged toward the Churchill- 
Chu line as the heat flux was increased. This slight shift to 
the left, although consistent during the three runs, was 
within the uncertainty in the region. Good repeatability was 
also found at the onset of nucleate boiling (indicated by 
arrows) and in the nucleate boiling region. Comparing the 
repeatability runs to the representative data set of Perry 
(Ref. 2], at low heat fluxes the data vary slightly but are 
within the uncertainty band. The data merge up to the point 
of incipience, but why Perry's run nucleated earlier than the 
present data is unclear since the runs were, in essence, 
identical. In the boiling region, Perry's data is slight wil mee 
the left of the present data but again lie within the 
uncertainty band for that region. 

In Figure 5.2, good repeatability is again shown POPE M 
heat fluxes. Comparing the three present repeatability runs 
to the representative data set of Perry [Ref. 2], again shows 
good reproducibility throughout. For heat fluxes < 1 kW/m’, 
the slopes of the data curves change, which indicates the 
transition from nucleate boiling back to natural convection. 
Based on Figures 5.1 and 5.2, data reproducibilst MM 
considered satisfactory with the single smooth tube work 


carried out by Perry (Ref. 2] on the same apparatus. 
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Figure 5.2 Decreasing Heat Flux Repeatability Runs using 
Smooth Tube Configuration CQ 


J? 


Figure 5.3 is a comparison of decreasing heat flux 
repeatability runs using an enhanced tube (TURBO-B) in 
configuration CO with the single tube data of Bertsch 
[Ref. 1]. For the present data, good repeatability is again 
obtained throughout. In comparing the current data with the 
data collected by Bertsch [Ref. 1], different boiling curves 
are easily seen. This was a result of two different TURBO-B 
test tubes being used. After examining both test tubes under 
a microscope (using 50x and 100x magnification), the TURBO-B 
test tube used in this study had wider notches between the 
'mushroom-like pedestals' (see Figure 3.4), indicating that a 
different fabrication process was used to make the tube. 
Above 20 kW/m^ on the plot, the current data curves are to the 
left of Bertsch's data (about 2 ?C less wall superheat), 
indicating better heat transfer. It is believed this is due 
to the wider notches in the tube surface allowing more bubbles 
to depart the tube surface at higher heat fluxes, removing 
more heat than the tube used by Bertsch. At about 10 kW/m’, 
the trend is reversed, the current data crosses and continues 
to fall below Bertsch's data (about 0.3 °C more wall 
Superheat). At low heat fluxes, the wider notches allow the 
refrigerant to flood back into the boiling channels (see 
Figure 3.4) and reduce the number of active sites, reducing 
the heat transfer from the tube. 

Figure 5.4 is a comparison of increasing heat flux 


repeatability runs using smooth tube configuration Cl with the 
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Figure 5.3 Decreasing Heat Flux Repeatability Runs using 
TURBO-B Tube Configuration CO 
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Figure 5.4 Increasing Heat Flux Repeatability Runs using 
Smooth Tube Configuration Cl with the Lower Tube Unheated 
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lower tube unheated and a single smooth tube data set 
(D0728ASCO) from Figure 5.1. The plot shows some scatter of 
@aeceat low heat fluxes that is attributed to the liquid pool 
mEuserature gradientS, as discussed previously. Again, good 
repeatability is found in the mid-range heat flux and boiling 
EE In comparison to the single tube data (configuration 
omand the Churchill and Chu correlation [Ref. 9], the upper 
tube heat transfer performance of the Cl configuration is 
clearly enhanced in natural convection. This enhancement is 
psebably due to the liquid pool surface being much closer to 
the upper test tube (approximately 10 mm) when it is position 
em increasing the circulation (i.e., local fluid velocity) 
around the upper tube and ultimately increasing heat transfer 
performance of the tube. Lake [Ref. 3], in his comparison of 
single tube R-114 data to his Cl configuration data, found the 
same enhancement. He concluded that the lower tube must be 
modifying the flow over the upper tube. In actuality, it is 
probably a combination of both the decreased distance to the 
pool surface and the lower tube modifying the flow over the 
upper tube that causes the enhancement. From this plot, it 
was concluded that the program DRPJY was performing properly. 

Figure 5.5 is the corresponding decreasing heat flux data 
runs for Figure 5.4. Good repeatability is again found for 
[I configuration runs. Comparing the single tube data 


u9728ASCO from Figure 5.2) to the Cl configuration data, 
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Figure 5.5 Decreasing Heat Flux Repeatability Runs using 
Smooth Tube Configuration Cl with the Lower Tube Unheated 
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indicates that in nucleate boiling with decreasing heat 
Mes, both tube configurations perform equally well. 

meguce sto 1s thé decreasing heat flux repeatability runs 
using TURBO-B tubes in configuration Cl with the lower tube 
unheated. Good repeatability is evident. Comparing the 
Single TURBO-B data set (DIOI2TBCO) from Figure 5.3 to the Cl 
data, at high heat fluxes both tube configurations perform 
equally well, as was for the smooth tube case (Figure 5.5). 
At low heat fluxes (< 8 kW/m’), the Cl configuration data is 
to the left of the CO configuration data, indicating an 
enhancement due to the addition of the lower tube and moving 
tetest tube closer to the surface of the liquid pool; this 
was not seen for the smooth tube case (Figure 5.5). As will 
be discussed later, the addition of a lower tube heat flux, 
shifts the Cl data curve even more towards the left. 

Based on Figures 5.1 through 5.6, the repeatability of the 
current data and reproducibility of previous single tube work 


performed by Perry [Ref. 2], is considered satisfactory. 


B. INFLUENCE OF THE LOWER TUBE 

To investigate the effects of a lower heated tube on an 
upper tube, five different tube configurations and fixed lower 
tube heat flux settings were used. The five tube 


EonbPIgurations selected were P/D ratios of 2.0, 1.8, 1.5, 1.2 
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Figure 5.6 Decreasing Heat Flux Repeatability Runs using 
TURBO-B Tube Configuration Cl with the Lower Tube Unheated 
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Eme 2 {rotated 30°). P/D of 2.0 was selected because it 
matches not only the vertical pitch used in naval shipboard 
flooded evaporators, but also the pitch u.ed in the flooded 
evaporator bundle apparatus at the Naval Postgraduate School 
Heat Transfer Laboratory. The fixed lower tube heat flux 
settings (LTHFS) used were 0, 1, 5, 10 and 20 kW/m’. In 
addition to these runs, data runs for each tube configuration 
were also done using a lower tube heat flux setting that was 
equal to the upper tube heat flux setting (UTHFS). Finally, 
Ara fixed UTHFS of 5, 10 and 20 kW/m’, the LTHFS was varied 
to investigate the influence of the strength of the convection 
effects from below. For each data run, the same tube surface 
mas used for the upper and lower tube; mixing of the tube 
types was not allowed. For the smooth tube runs, the upper 
tube was started in the natural convection region at low heat 
fluxes, was run up and down through various heat flux 
Settings, and was returned to the natural convection region. 
For the TURBO-B tubes, only decreasing upper tube heat flux 
runs were recorded to prevent overlapping/repeating data (the 
TURBO-B tubes did not enter the natural convection region at 
low heat fluxes). 

To evaluate the effect of a lower heated tube on an upper 
tube, the data is presented by the effect of changing the P/D 
ratio for a common LTHFS, the effect of varying the LTHFS for 
a common P/D ratio and the effect of varying the lower tube 


heat flux for a fixed upper tube heat flux setting. 
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l. The Effect of Tube Spacing 
a. Smooth Tubes with Lower Tube Unheated 

Figure 5.7 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with the lower 
tube unheated and the corresponding Churchill and Chu 
correlation [Ref. 9]. As seen in the repeatability runs, at 
low heat fluxes the data is shifted to the left of the 
Churchill-Chu line and as the heat flux is increased, the data 
tends to approach the Churchill-Chu line up to the incipient 
boiling point (indicated by arrows). This implies that the 
Circulation occurring within the evaporator enhances the 
natural convection that occurs from a single tube immersed in 
an infinite expanse of fluid. Following nucleation of the 
upper tube, the boiling curve for P/D of 1.8 is to the left of 
the others and appears to be yielding the best heat transfer 
performance. It is unclear why this is so, because it is so 
repeatable and the same smooth test tubes were used. Lake 
[Ref. 3] also found that P/D of 1.8 gave the best heat 
transfer performance. 

Figure 5.8 shows the corresponding decreasing heat 
flux data for Figure 5.7. With decreasing heat flux, a P/D of 
1.8 again provides the best heat transfer performance until 


the lowest heat flux region, as did Lake's data [Ref. 3]. 
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Figure 5.7 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with the Lower Tube Unheated 
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Figure 5.8 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with the Lower Tube Unheated 
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b. Smooth Tubes with LTHFS of 1 kW/n 

Hicure Shows Ene) jncreasing heat flux data 
taken for the five tube configurations (Cl-C5) with a heat 
flux of 1 kW/m’ on the lower tube. Both the upper and lower 
tubes for all configurations started in the natural convection 
region. At low heat fluxes, some scatter is present which 
becomes smaller as the data approaches the incipient boiling 
Pommes and tends towards the Churchill and Chu correlation 
[Ref. 9], as seen previously. Note that the data for a P/D of 
1.2 crosses the Churchill-Chu line and appears to nucleate 
earlier than the others. For this small spacing, the heated 
plume from the lower tube can actually reduce natural 
convection from the top tube since the local ambient pool 
temperature is increased. This creates a larger wall 
temperature anda larger wall superheat (Twall - Tsat). After 
nucleation, a P/D of 1.8 again provides better heat transfer 
in the boiling region than the other P/D ratios. This trend 
Suggests that in the boiling region the enhancement provided 
by a P/D of 1.8 is due to the tube spacing itself and not from 
heating the lower tube. 

Figure 5.10 is the decreasing heat flux data for to 
Figure 5.9. A P/D of 1.8 provides the best heat transfer 
performance until approximately 4 kW/m^, which is where the 
boiling curves seem to merge. With a LTHFS of 1 kW/m’, it 


appears to force the decreasing boiling curves to merge at a 
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Figure 5.9 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with 1 kW/m* on the Lower Tube 
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Figure 5.10 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with 1 kW/m” on the Lower Tube 
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higher heat flux (4 kW/m^) than with the lower tube unheated 
(Sa 
c. Smooth Tubes with LTHFS of 5 kW/m 

Figure 5.11 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with 5 kW/m’ on 
the lower tube. Both upper and lower tubes for all the 
configurations started in the natural convection region. At 
low heat fluxes, the data is to the left of the Churchill and 
Chu correlation [Ref. 9], as seen in Figures 5.7 and 5.9. As 
the heat flux is increased, some interesting things happen. 
The data for a P/D of 1.2 {rotated 30º) quickly cross sP 
Churchill-Chu line and begins to parallel it, then nucleates 
at a low heat flux of about 3 kW/m’. This is believed to be 
an anomaly and not an effect. Of the remaining P/D ratios, 
1.8 is the only one to tend toward the Churchill-Chu linemirn 
increasing heat flux. The P/D ratios of 2.0, 1.5 and 1.2 seem 
to parallel the Churchill-Chu line, probably due to the L-HFS 
intensifying the forced convection around the upper tube. In 
all cases, the lower tube remained in natural convection until 
the upper tube nucleated. Once the upper tube nucleated, the 
lower tube also nucleated. In the boiling region, a P/D of 
1.8 again gave the best enhancement. 

Figure 5.12 is the decreasing heat flux data for 


Figure 5.11. Above 10 kW/m’, P/D of 1.8 is clearly providing 
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Figure 5.11 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with 5 kW/m^ on the Lower Tube 
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Figure 5.12 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with 5 kW/m* on the Lower Tube 
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the best heat transfer. Below 10 kW/m’, the upper tube is 
returning to natural convection, but the lower tube is still 
nucleating. The bubbles departing the lower tube impinge upon 
the upper tube causing a significant increase in heat transfer 
performance, from bubble sweeping, compared to the data for a 
LTHES of 0 and 1 kW/mí. Also, a LTHFS of 5 kW/m? eliminates 
any upper tube hysteresis effect. 
d. Smooth Tubes with LTHFS of 10 kW/m 

Figure 5.13 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with 10 kW/m’ 
on the lower tube. For this LTHFS, the lower tube nucleated 
at the start of the run for a P/D of 1.5, but the upper tube 
remained in natural convection. The effect of the bubbles 
from the lower tube impacting and sweeping over the upper tube 
is seen to significantly enhance the heat transfer 
performance of the upper tube in the low heat flux regions. 
For the other P/D ratios, the incipient boiling points were 
very random; at the largest spacing, nucleation occurred 
quickly, whereas at the smallest spacing, nucleation was 
delayed. Possibly, the largest spacing (P/D of 2.0) allowed 
the heated plume from the lower tube to thin the boundary 
layer around the upper tube, but the closest spacing (P/D of 
1.2) enabled the heated plume to 'skirt' around the upper tube 
and not disturb the upper tube boundary layer as much, 


delaying the point of incipience. Again, once the upper tube 
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Figure 5.13 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with 10 kW/m* on the Lower Tube 
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nucleated, the lower tube also nucleated shortly thereafter. 
In the fully developed nucleate boiling region, a P/D of 1.8 
still gives the best performance. 

Figure 5.14 is the decreasing heat flux data for 
Figure 5.13. A P/D of 1.8 provided the best heat transfer 
until the heat flux was lowered to about 15 kW/m’, where no 
gain over the other P/D ratios is noticed. This concurs with 
the findings of Lake [Ref. 3]. Below 10 kW/m”, the nucleating 
lower tube is definitely enhancing the upper tube performance 
mou! P/D ratios, although, the P/D of 1.2 (rotated 309) is 
slightly to the right of the other curves. This is probably 
caused by the angular offset of the upper tube, not allowing 
the bubbles from the lower tube to directly impinge upon the 
upper tube. Therefore, a 30° rotation slightly decreases the 
performance of the upper tube compared to an in-line tube 
configuration. 

e. Smooth Tubes with LTHES of 20 kW/m 

Figure 5.15 shows the increasing heat flux data 
taken for the five tube configurations (C1-C5) with 20 kW/m' 
on the lower tube. In each run, the lower tube was nucleating 
at the beginning of the run and remained nucleating 
throughout. The only effect from P/D spacing is noticed above 
20 kW/m? where a P/D of 1.8 again provides the best 


performance. 
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Figure 5.14 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with 10 kW/m* on the Lower Tube 
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Figure 5.15 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with 20 kW/m? on the Lower Tube 
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Figure 5.16 is the decreasing heat flux data for 
Figure 5.15. Again, below 20 kW/m?, only a P/D NA 
(rotated 30?) drops slightly below the other curves at mid- 
range heat fluxes. It appears that this LTHFS enables more 
bubbles from the lower tube to impinge upon the upper tube, 
compared to a LTHFS of 10 kW/m’. 

f. Smooth Tubes with UTHFS equal LTHFS 

Figure 5.17 shows the increasing heat flux data 
taken for the five tube configurations (Cl-C5) with the UTHFS 
equal to the LTHFS throughout the entire length of ChE moe 
At the beginning of all the runs, both tubes were in the 
natural convection region. In the natural convection region, 
a P/D of 1.2 (rotated 30°) is the only configuration PME 
very close to the Churchill and Chu correlation [R° mE IE 
This was not expected nor can be explained. Again, a random 
incipient point was obtained and in the boiling region, a P/D 
of 1.8 was the best performer. 

Figure 5.18 is the decreasing heat flux data for 
Figure 5.17. Above 20 kW/m’, a P/D of 1.8 still provides the 
best performance and below 20 kW/m? it appears that as the P/D 
decreases, so does the enhancement of the upper tube. Between 
9 and 5 kW/m’, it's believed that the tubes shift fromm@ene 
boiling region back to natural convection. This may causem 


decrease of enhancement because, as speculated earlier, as the 
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Figure 5.16 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux VEM 20 kW/m* on the Lower Tube 
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Figure 5.17 Comparison of Smooth Tube P/D Ratios for 
Increasing Heat Flux with UTHFS equal to LTHFS 
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Figure 5.18 Comparison of Smooth Tube P/D Ratios for 
Decreasing Heat Flux with UTHFS equal to LTHFS 
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P/D is decreased so too are the convection effects from the 
lower meUDe: 
g. TURBO-B Tubes with Lower Tube Unheated 

Figure 5.19 shows the decreasing heat flux data 
taken for the five tube configurations (Cl-C5) with lower tube 
unheated. At the start of each run, an undetermined amount of 
nucleation was occurring from the lower tube due to the 
vigorous boiling action off the upper tube at the run start 
heat flux (100 kW/m’). A very consistent boiling curve above 
10 kW/m^ indicates that there is no enhancement from the lower 
tube at any P/D ratio. Below 5 kW/m^, the rotated 
configuration data agrees well with the single tube data in 
Page 559. This indicates that the lower tube is not 
modifying the flow around the upper tube and the upper tube is 
performing as a single tube which provides less heat transfer 
than the in-line configurations. 

h. TURBO-B Tubes with LTHFS of 1 kW/m’ 

Figure 5.20 shows the decreasing heat flux data 
taken for the five tube configurations (C1-C5) with 1 kW/m’ on 
the lower tube. At the start of each run, both the upper and 
lower tubes were nucleating. Very little change is noticed 
between Figures 5.19 and 5.20. Above 20 kW/m’, the plots are 
identical. Below 10 kW/m’, the LTHFS of 1 kW/m? moved the data 


curves slightly to the left. This indicates that the greater 
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Figure 5.19 Comparison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with the Lower Tube Unheated 
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Figure 5.20 Comparison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with 1 kW/m* on the Lower Tube 
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the LTHFS, the more enhanced the mid-range and lower heat flux 
Medion becomes. As seen in the previous figure, a P/D of 1.2 
(rotated 30?) is again less enhanced than the other P/D 
ratios. 

i. TURBO-B Tubes with LTHFS of 5 kW/n’ 

Figure 5.21 shows the decreasing heat flux data 
taken for the five tube configurations (C1-C5) with 5 kW/m’ on 
the lower tube. At the start of each run, both tubes were 
vigorously boiling. Above 20 kW/m’, the data is consistent 
with the two previous plots, but below 10 kW/m’, a definite 
separation between a P/D of 1.2 (rotated 30?) and the other 
P/D ratios is seen. Also for heat fluxes < 10 kW/m’, the in- 
line configurations provide better heat transfer than the same 
Configurations seen in Figure 5.20. This indicates that as 
LTHFS increases, more bubbles depart the lower tube and 
increase the upper tube heat transfer performance for the in- 
line configurations. or the rotated configuration, the 
bubbles must be departing the lower TURBO-B tube differently 
(compare to a lower smooth tube) and are not impinging upon 
the upper TURBO-B tube. This explains the approximately 40 - 
45% decrease in heat transfer performance for the rotated 
eom guration compared to the in-line configurations. 

j. TURBO-B Tubes with LTHFS of 10 kW/m 
Figure 5.22 shows the decreasing heat flux data 


taken for the five tube configurations (C1-C5) with 10 kW/m‘ 
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Figure 5.21 Comparison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with 5 kW/m^ on the Lower Tube 
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Figure 5.22 Comparison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with 10 kW/m^ on the Lower Tube 
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on the lower tube. Again, both tubes were nucleating at the 
beginning of each run. There are no obvious changes above 20 
kW/m? from the previous plots. Below 10 kW/m’, all the curves 
seem to have shifted to the left compared to the previous 
case. 
k. TURBO-B Tubes with LTHFS of 20 kW/né 
Figure 5.23 shows the decreasing heat flux data 
taken for the five tube configurations (C1-C5) with 20 kW/m’ 
on the lower tube. Figures 5.22 and 5.23 show simim, 
results, therefore, increasing the LTHFS from 10 to 20 kW/m’ 
provides a small change in upper tube heat transfer 
performance. 
l. TURBO-B Tubes with UTHFS equal LTHFS 
Figure 5.24 shows the decreasing heat finu 
taken for the five tube configurations (C1-C5) with the UTHFS 
equal to the LTHFS. There is not much change in this plot 
relative to the trend set by the previous plots. 
2. The Effect of a Fixed Lower Tube Heat Flux 
a. Smooth Tubes with P/D of 2.0 
Figure 5.25 shows the increasing heat flux data 
taken for a P/D of 2.0 at five different heat flux settings on 
the lower tube. The Churchill and Chu correlation [Ref. 9] is 
also shown for comparison. For a LTHFS of 20 kW/m’, the lower 
tube had nucleated immediately at the beginning of the run and 


remained nucleating for the length of the run. All other 
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Figure 5.23 Comparison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with 20 kW/m? on the Lower Tube 
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Figure 5.24 Comparison of TURBO-B Tube P/D Ratios for 
Decreasing Heat Flux with the UTHFS equal to the LTHFS 
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Figure 5.25 Comparison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Tubes with P/D of 2.0 
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LTHFS data runs, started in natural convection. Ine 
natural convection region, the data are shifted to the left of 
the Churchill and Chu correlation [Ref. 9], but approach" 
heat flux is increased. This is probably caused Oy ES m 
circulation patterns within the small evaporator enhancing the 
heat transfer mentioned previously. In natural convection, 
the effect of a LTHFS had little influence on the upper tube 
heat transfer behavior; in fact, this is repeated for all the 
P/D ratios. The point of incipience was very repeatable, with 
the exception of a LTHFS of 10 kW/m’. A heat flux of 10 kW/m? 
is approximately the heat flux that 'trips' the nucleation 
sites causing boiling to occur from a smooth tube, so shortly 
after the run started, the lower tube nucleated and caused the 
upper tube to nucleate quickly. At high heat fluxes, all the 
data falls on one nucleate boiling curve indicating that the 
LTHFS has little effect on the heat transfer from the upper 
tube in this region. 

Figure 5.26 shows the decreasing heat flux data for 
the data run for Figure 5.25. Important to note that the 
upper tube starts and finishes each run in natural convection, 
but the lower tube might not. For a LTHFS of 0 and 1 kW/m', 
the lower tube starts and finishes in natural convection; for 
a LTHFS of 20 kW/m’, it starts and finishes in nucleate 
Doiling., and for a LTHFS of 5 and 10 kW/m’, it star NEN 
natural convection and finishes in nucleate boiling. At high 


heat fluxes, all the data lay on one nucleate boiling curve as 
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Figure 5.26 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Smooth Tubes with P/D of 2.0 
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expected. Below 20 kW/m’, five separate curves are seen. The 
curves with a LTHFS of O and 1 kW/m? return to the natural 
convection region, indicating a deactivation of the nucleation 
Sites as the upper tube heat flux is decreased, showing the 
hysteresis effect of the upper tube. At a heat flux below 5 
kW/m’, the curves with a LTHFS of 5, 10 and 20 kW/m’ begin to 
parallel the natural convection curve found on the increasing 
heat flux plot and have significant enhancement over the O and 
1 kW/m^ curves at low heat fluxes. Also, the hysteresis 
effect is eliminated by the vigorously nucleating lower tube. 
This same phenomena was found by Lake [Ref. 3]. He attributed 
this enhancement to strong convection effects on the upper 
tube from the bubbles nucleating off the lower tube. 
Therefore, significant enhancement of heat transfer of an 
upper tube is possible when the lower tube is nucleating and 
very little enhancement is noticed when both tubes are in 
natural convection. 
b. Smooth Tubes with P/D of 1.8 

Figure 5.27 shows the increasing heat flux data 
taken for a P/D of 1.8 at five lower tube heat flux settings. 
For a LTHFS of 20 kW/m’, the lower tube had nucleated 
immediately at the beginning of the run and remained 
nucleating for the entire length of the run. All other han: 
started in natural convection. In the natural convection 


region, the data are shifted to the left of the Churchill and 
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Figure 5.27 Comparison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Tubes with P/D of 1.8 
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Chu correlation [Ref. 9], but approach it as heat MEN. 
increased, as was the case with a P/D of 2.0. Again, the 
effect of the LTHFS had little influence of the upper tube 
heat transfer behavior. The point of incipience was again, 
very repeatable. This time, unlike for a P/D of 2.0, the 
lower tube remained in natural convection until the upper tube 
nucleated. Again, at high heat fluxes, all the data falls on 
one nucleate boiling curve indicating that the LTHFS has 
little effect upon the heat transfer in this region. 

Figure 5.28 shows the decreasing heat flux data for 
the data runs for Figure 5.27. The data trend is similar to 
that found in Figure 5.26. Therefore, significant enhancement 
of heat transfer from an upper tube is again seen when the 
lower tube is nucleating with larger than 5 kW/m? on the lower 
tube. 

c. Smooth Tubes with P/D of 1.5 

Figure 5.29 shows the increasing heat flux data 
taken for a P/D of 1.5 at five lower tube heat flux settings. 
For a LTHFS of 10 and 20 kW/m’, the lower tube had nucleated 
immediately at the beginning of the run and remained 
nucleating for the length of the run. All other LTHES started 
in natural convection. This plot shows similar results ec 


on Figure 5.27 and most of the comments still apply. 
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Figure 5.28 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Smooth Tubes with P/D of 1.8 
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Figure 5.29 Comparison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Tubes with P/D of 1.5 
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Figure 5.30 shows the decreasing heat flux data for 
the data runs shown in Figure 5.29. Again, similar results 
are obtained as before. Significant enhancement in the heat 
transfer from an upper tube is obtained when the lower tube is 
nucleating with a heat flux larger than 5 kW/m’. 

d. Smooth Tubes with P/D of 1.2 

Figure 5.31 shows the increasing heat flux data 
taken for a P/D of 1.2 at five lower tube heat flux settings. 
Homma’ LTHFS of 20 kW/m’, the lower tube had nucleated 
immediately at the beginning of the run and remained 
nucleating for the length of the run. All other LTHFS started 
in natural convection. This plot shows similar results seen 
on Figure 5.29 and most of the comments still apply. 

Figure 5.32 shows the decreasing heat flux data for 
the data runs shown in Figure 5.31. This plot shows similar 
results seen on Figure 5.30 and again, most of the comments 
still apply. 

e. Smooth Tubes with P/D of 1.2 (Rotated 30°) 

Figure 5.33 shows the increasing heat flux data 
taken for a P/D of 1.2 (rotated 30°), at five lower tube heat 
flux settings. For a LTHFS of 20 kW/m”, the lower tube had 
nucleated immediately at the beginning of the run and remained 
musEcatrng for the length of the run. All other LTHFS started 
in natural convection. The increasing data for this 


configuration is more concentrated than the others seen 
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Figure 5.30 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Smooth Tubes with P/D of 1.5 
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Figure 5.31 Comparison of Lower Tube Heat Flux Settings on 
Increasing Heat Flux of Smooth Tubes with P/D of 1.2 
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Figure 5.32 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Smooth Tubes with P/D of 1.2 
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Figure 5.33 Comparison of Lower Tube Heat Flux Settings on 
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Increasing Heat Flux of Smooth Tubes with P/D of 1.2 


(Rotated 30°) 


previously. This might be due to the offset orientar EE i 
the upper tube causing the data to be less scattered. 

Figure 5.34 shows the decreasing heat flux data for 
the data runs shown in Figure 5.33. At high heat fluxes, all 
the data lay on one nucleate boiling curve. Below 10 kW/m’, 
separate curves are formed and as Seen before. Significant 
enhancement of heat transfer from an upper tube is obtained 
when the lower tube is nucleating with a heat flux larger than 
5 kW/m’. 

f. TURBO-B Tubes with P/D of 2.0 

Figure 5.35 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
At the start of each run, a very consistent boiling curve is 
seen above 20 kW/m’, indicating that above a UTHFS of 20 kW/m’, 
there is no enhancement from increasing the LTHFS. Below 10 
kW/m’, raising the LTHFS does increase the heat transfer 
performance of the upper tube. At an upper tube heat flux of 
2 kW/m’, an increase in performance of 2.4 is found going from 
an unheated lower tube to a LTHFS of 20 kW/m’, compared to an 
increase of 3.6 for the smooth tube array (Figure 5.26). 

g. TURBO-B Tubes with P/D of 1.8 

Figure 5.36 shows the decreasing heat flux data for 

the TURBO-B tubes taken at five lower tube heat flux settings. 


Similar results are found compared to the P/D of 2.0 case. At 
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Figure 3.34 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of Smooth Tubes with P/D of 1.2 
(Rotated 309?) 
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Figure 5.35 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 2.0 
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Figure 5.36 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 1.8 
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an upper tube heat flux of 2 kW/m*, an increase in performance 
of 1.6 is found going from a unheated lower tube to a LTHFS of 
20 kW/m*, compared to an increase of 3.7 for the smooth tube 
array (Figure 5.28). Unlike the smooth tube case where a P/D 
of 1.8 produced the best heat transfer performance in nucleate 
boiling, no optimal P/D ratio was found for the TURBO-B tube 
array. 
h. TURBO-B Tubes with P/D of 1.5 

Figure 5.37 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
Similar results are found compared to the P/D of 1.8 case. At 
an upper tube heat flux of 2 kW/m^, an increase in performance 
of 2.2 is found going from a unheated lower tube to a LTHFS of 
20 kW/m*, compared to an increase of 3.5 for the smooth tube 
array (Figure 5.30). 

i. TURBO-B Tubes with P/D of 1.2 

Figure 5.38 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
Similar results are found compared to the P/D of 1.5 case. At 
an upper tube heat flux of 2 kW/m*, an increase in performance 
of 2.2 is found going from a unheated lower tube to a LTHFS of 
20 kW/m^, compared to an increase of 5.4 for the smooth tube 


array (Figure 5.32). 


106 








< 

> 4 

E 10 

= 

x= 

> 

OG 

e 

D 

o 

o 

La 

3 

= 

a. 10? 

= 

- | 
lower tube 0 kW/m'2 
lower tube 1 kW/m'2 
lower tube 5 kW/m'2 
lower tube 10 kW/m' 2 
lower tube 20 kW[m'2 
10? 


107 10° 10! 
Wall Superheat/(K) 


Figure 5.37 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 1.5 
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Figure 5.38 Comparison of Lower Tube Heat Flux Settings on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 1.2 
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j. TURBO-B Tubes with P/D of 1.2 (Rotated 30?) 

Figure 5.39 shows the decreasing heat flux data for 
the TURBO-B tubes taken at five lower tube heat flux settings. 
A very consistent boiling curve is again seen above 10 kW/m’. 
Below 10 kW/m’, the data seems less scattered and slightly to 
BH Ncrt of any other configuration, indicating that rotating 
the upper tube reduces the enhancement produced from 
increasing the LTHFS in the low heat flux regions. At an 
upper tube heat flux of 2 kW/m’, an increase in performance of 
1.7 is found going from a unheated lower tube to a LTHFS of 20 
kW/m’, compared to an increase of 3.5 for the smooth tube 
array (Figure 5.34). 

3. The Effect of a Variable Lower Tube Heat Flux 
a. Smooth Tubes with P/D of 2.0 

Figure 5.40 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is varied. The three fixed UTHFS used were 5, 10 
and 20 kW/m’. The labeling in the plot legend describes which 
data were taken uSing increasing or decreasing lower tube heat 
flux. The sharp transitions in the increasing curves are the 
points where the lower tube nucleated and caused the upper 
tube to shift from natural convection to nucleate boiling 
(visually verified during the data run), significantly 
increasing the upper tube heat transfer coefficient, as 


expected. The maximum upper tube heat transfer coefficient 
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Figure 5.39 Comparison of Lower Tube Heat Flux Setting on 
Decreasing Heat Flux of TURBO-B Tubes with P/D of 1.2 
(Rotated 30º) 
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Figure 5.40 Heat Transfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 2.0 
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for all three UTHFS was 1850 W/m’:K, which indicates ames 
a lower tube heat flux of 80 kW/m’, the heat transfer 
performance of the upper tube is independent of the lower tube 
heat flux. The decreasing curves show that for an UTHFS of 20 
kW/m’, the lower tube heat flux has little influence on the 
upper tube heat transfer. For the UTHFS of 10 kW/m’, as the 
lower tube heat flux is decreased, the upper tube remains in 
nucleate boiling, but for the UTHFS of 5 kW/m’, the upper tube 
returns to the natural convection region as the lower tube 
heat flux is decreased. 
b. Smooth Tubes with P/D of 1.8 

Figure 5.41 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is varied. Similar results were obtained as in the 
P/D of 2.0 case, except that the curves did not merge at the 
highest heat flux setting. The maximum heat transfer 
coefficient for the UTHFS of 5 and 10 kW/m? converged at 
approximately 1700 W/m’:K and for the UTHFS of 20 kW/m’ the 
value was about 1900 W/m*:*K. It is not known why the curves 
did not merge at high lower tube heat fluxed. Comparing these 
results to the P/D of 2.0 results, at the smallest lower tube 


heat fluxes, the P/D of 1.8 provides better heat transfer. 
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Figure 5.41 Heat Transfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 1.8 
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Figure 5.42 Heat Transfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 1.5 
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c. Smooth Tubes with P/D of 1.5 
Figure 5.42 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is varied. The maximum heat transfer coefficient 
for all three of the UTHFS converged at approximately 1800 
W/m^:K, similar to the P/D of 2.0 case. The decreasing curves 
show very similar behavior found at a P/D of 2.0. Comparing 
these results to data of a P/D of 2.0 and 1.8, a P/D of 2.0 
and 1.5 are similar in the amount of enhancement achieved, 
whereas a P/D of 1.8 outperforms them both, in the mid to low 
heat flux regions. 
d. Smooth Tubes with P/D of 1.2 
Figure 5.43 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is increased. The maximum heat transfer coefficient 
for all three of the UTHFS converged, as did the P/D of 2.0 
and 1.5, at approximately 1990 W/m^:K, highest value of all 
the configurations tested. Therefore, the best heat transfer 
performance was given by a P/D ratio of 1.8 when the lower 
tube heat flux was held constant, but a P/D ratio of 1.2 when 
the upper tube heat flux was held constant. No decreasing 
curves were recorded for these data runs. 
e. Smooth Tubes with P/D of 1.2 (Rotated 30?) 
Figure 5.44 shows the heat transfer performance of 


an upper tube, held at a constant heat flux, as the lower tube 
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Figure 5.43 Heat Transfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 1.2 
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Figure 5.44 Heat Transfer Performance of a Smooth Upper 
Tube at a Constant Heat Flux for a P/D of 1.2 (Rotated 30º) 


ie 


heat flux is increased. As seen for the P/D of 1.8 case, the 
curves do not merge at the highest lower tube heat flux. The 
maximum heat transfer coefficient for the two lower curves 
converged at about 1500 W/m^:K and the upper curve was 1725 
W/m^-K. Comparing Figures 5.43 and 5.44, for an UTHE SMA 
kW/m?, the data curves are identical for lower tube heat 
fluxes up to 8 kW/m’, indicating that the forced conve. Coua 
effects from the lower tube have very little influence on the 
upper tube heat transfer behavior. Above lower tube heat 
fluxes of 10 kW/m*, the in-line configuration provides better 
heat transfer performance because the lower tube is in 
nucleate boiling and all of the bubbles are impacting the 
upper tube providing enhancement by bubble pumping; as 
opposed to the offset configuration, where most of the bubbles 
coming off the lower tube miss the upper tube. The same idea 
holds true for the other data curves. The  in-line 
configurations are more enhanced due to bubble pumping from 
the lower tube. No decreasing curves were recorded for these 
data runs. 
f. TURBO-B Tubes with P/D of 2.0 

Figure 5.45 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. At high lower tube heat fluxes, all 
three of the UTHFS appear to perform almost equally. For an 


UTHFS of 20 kW/m’, the lower tube heat flux has little 
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Figure 5.45 Beat Transfer Performance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 2.0 
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influence on the upper tube heat transfer and for an UTHES TOT 
5 and 10 kW/m”, as the lower tube heat flux is decreased so 
follows the performance of the upper tube (more for 5 than 10 
kW/m‘) . 
g. TURBO-B Tubes with P/D of 1.8 
Figure 5.46 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. Comparing these curves to the ones 
for a P/D of 2.0, a small enhancement in performance is gained 
in the mid and low heat flux regions using a P/D of 1.8, which 
is an expected result. 
h. TURBO-B Tubes with P/D of 1.5 
Figure 5.47 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. There is very little difference in 
performance between a P/D of 1.5 and 1.8. Comparing the P/D 
of 1.5 to the 2.0 data, a small decrease in performance is 
noted for lower tube heat fluxes > 20 kW/m’ using the P/D of 
1.5 configuration, which cannot be explained. 
i. TURBO-B Tubes with P/D of 1.2 
Figure 5.48 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. Similar performance is noted for a 
P/D ratio of 1.2 and 1.5, therefore, the comments made 


previously still apply. 
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Figure 5.46 Heat Transfer Performance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 1.8 
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Figure 5.47 Heat Transfer Performance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 1.5 
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Figure 5.48 Heat Transfer Performance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 1.2 
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j. TURBO-B Tubes with P/D of 1.2 (Rotated 30?) 

Figure 5.49 shows the heat transfer performance of 
an upper tube, held at a constant heat flux, as the lower tube 
heat flux is decreased. Comparing Figures 5.48 and 5.49, for 
an UTHFS of 20 kW/m’, the data curves are almost identical, 
indicating that the lower tube orientation and heat flux have 
very little influence on the upper tube heat transfer 
behavior. For the other two data curves, the in-line 
configuration provides better heat transfer performance 
because of bubble pumping, discussed earlier. For the offset 
configuration, all of the data curves do not merge and are 
almost horizontal lines, indicating there is no enhancement by 
bubble pumping. It is believed that most of the bubbles 
coming off the lower tube miss the upper tube because the 
bubbles depart the TURBO-B tube surface at the upper most 
section of the tube, which is different than how bubbles 
depart the smooth tube surface. This phenomena needs to be 
further investigated. No decreasing curves were recorded for 


these data runs. 


C. R-124/R-114 PERFORMANCE COMPARISON 

Figure 5.50 compares the decreasing heat flux data for R- 
124 found in this study with the R-114 data of Lake [Ref. 3] 
for smooth tube configurations listed in the plot legend, with 


the lower tube unheated. In the mid to high heat flux 


124 
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Lower Tube Heat Flux/(W/m 2) 


Figure 5.49 Heat Transfer Performance of a TURBO-B Upper 
Tube at a Constant Heat Flux for a P/D of 1.2 (Rotated 30 ) 
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regions, the P/D ratios of 1.8 and 1.5 demonstrated improved 
performance with the R-124, but with a P/D of 2.0, little or 
no improvement was found. In the low heat flux region, the 
uncertainty is too large to make any definitive conclusions. 

Figure 5.51 compares the decreasing heat flux data for R- 
124 with the R-114 data of Lake [Ref. 3] for the smooth tube 
configurations listed in the plot legend, with 10 kW/m’ on the 
lower tube. In the high heat flux region, a P/D of 1.8 and 
1.5 again provided better performance using R-124, but for R- 
114, a P/D of 2.0 displayed slightly better performance. In 
the  mid-range heat fluxes, R-114 exhibited improved 
performance in each of the three configurations. Again, in 
the low heat flux region, there was too much uncertainty in 


the data to draw any viable conclusions. 
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VI; CONCLUSIONS 


The interaction effects of a lower heated tube on pool 
boiling of R-124 from an upper horizontal tube has been 
studied for tube two types: smooth and TURBO-B. Based on the 


data obtained, the following conclusions are made: 


1. In natural convection, the effect of a lower heated tube 
on the heat transfer from an upper tube was small. 


2. For a smooth tube array, a vigorously nucleating lower 
tube eliminated upper tube hysteresis. 


3. For a smooth tube array in nucleate boiling, P/D of 1.8 
gave the best upper tube heat transfer performance. NO 
optimal P/D was found for the TURBO-B tube array. 


4. For both arrays with upper tube heat fluxes « 5 kW/m, 
the enhancing effect of bubbles from a lower tube was 
dramatic. This enhancement increased as lower tube heat 
flux increased. When upper tube heat fluxes were > 20 
kW/m’, all enhancement disappeared. 


5. For a TURBO-B tube array, a 30 degree offset of the upper 
tube reduced the upper tube heat transfer performance 
(Compared to the in-line configurations). This may 
indicate bubbles depart TURBO-B tubes differently than 
smooth tubes. 


6. For a smooth tube array with the lower tube unheated and 
an upper tube heat flux of > 3 kW/m’, the performance of 
R-124 compared to R-114, improved as the P/D ratio was 
reduced. With a nucleating lower tube (at 10 kW/m’), 
again the performance of R-124 improved as the P/D ratio 
decreased, but only for upper tube heat fluxes of > 40 
kW/m‘. 
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VII. RECOMMENDATIONS 


Conduct further studies using TURBO-B tubes to verify 
that bubbles do depart the TURBO-B tubes differently than 
smooth Lubes. 


Conduct similar tests using different enhanced tubes to 
better understand the influence of a lower heated 
enhanced tube on the heat transfer behavior of an upper 
tube. 


. Modify the existing evaporator to accept a larger number 


of tubes (three and four smooth or enhanced tubes) and 
conduct tests to study the effect of multiple heated 
lower tubes on the heat transfer performance of an upper 
tube. 


Conduct similar tests using multiple tube configurations 


with R-124/oil mixtures to determine the effect of oil on 
the heat transfer performance of an upper tube. 


A 


APPENDIX A. PROGRAM SETUPJY 


Program SETUPJY was used to evaluate the system 


performance prior to commencing a data collection run. This 


program is written in Hewlett-Packard Basic 5.0 language and 


enables the operator to: 


Co N 


Monitor the coolant sump temperature. 

Monitor the evaporator liquid pool average temperature. 
Monitor all the thermocouple channel output temperatures. 
Monitor the voltage, current and power supplied to the 
upper tube heater (main heater) as well as the heater 
element electrical resistance. 

Monitor the voltage, current and power supplied to the 


lower tube heater (auxiliary heater) as well as the 
heater element electrical resistance. 


lo] 


10! 
20! 
301 
40! 
50! 
60! 
70 

80 


90 
100 
110 
120 
150 
160 
170 
180 
190 
200 
210 
230 
240 
250 
260 
270 
280 
290 
300 
310 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
5101! 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 


PROGRAM: SETUPJY 


DATE: 


MAY 15, 
LT DEAN SUGIYAMA 


PROGRAMMER: 


MODIFIED 3Y LANNIE LAKE JAN 22, 
MODIFIED BY GEORGE D. PERRY 10 APR 
MODIFIED 3Y JOE YUSICIAN MAY 15, 


1993 


1997 
1993 
1993 


COM /Cc/ C(7) 

DATA 0.10086091,25727.94369,-767345.8295,78025595.81,-9247486589,6.97688E+ 
11,-2.66192E*13 
DATA 3.94078E*14 
READ C(*) 
ON XEY 1,15 GOTO 150 


2RINTER IS 1 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
BEEP 

INPUT 


USING 
USING 
USING 
USING 
USING 
USING 
USING 


Ido 


"LX, “SELECT OPTION” ** 
"6X,""QzMONITOR SUMP""" 

"6X,"" 13MONITOCR LIQUID""" 
"6X,""23CHBCK THERMCCOUPLES""" 
"6X,""3zCHECK MAIN HEATER""" 
"6X,""4=CHECK AUX HEATERS""" 
"6X,""5-EXIT PROGRAM""" 


IF 
IF 
IF 
IF 
IP 
LE 
IF 


Ido»5 
Ido=0 
Idos1 
Ido=2 
Ido=3 
Ido=4 
Idos5 


PRINT 
PRINT "SUMP TEMPERATURE DEG C " 
PRINT 
OUTPUT 709;"AR AF19 AL19 VRS" 
OUTPUT 709;"AS SA" 

SumsQ 
FOR Jz1 TO 5 
ENTER 709;E 
SumsSum*& 
NEXT J 
Bave=Sum/5 

Tamp=FNTvsv( Save) 

PRINT USING "4X,MDD.DD";Temp 
BEEP 
PRINT 
WAIT 5 
GOTO 340 


PRINT 
PRINT "LIQUID TEMPERATURE DEG C” 
PRINT 
OUTPUT 709;"AR AF16 AL18 VR5” 
Sum=0 
FOR Is1 TO 2 
OUTPUT 709;"AS SA” 
ENTER 709;B 
Sum=Sum+E 

NEXT I 
BaveeSum/2 

TempeFNTvsv(Eave) 

PRINT USING "4X,MDD.2D";Temp 
BEEP 
WAIT 5 


GOTO 610 


Idos5 
340 
520 
690 
870 
870 
1170 


THEN 
THEN 
THEN 
THEN 
THEN 
THEN 
THEN 


I 


580! 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860! 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 


PRINT 

PRINT "TEMPERATURE DEG C" 
OUTPUT 709;"AR AF00 AL19 VRS" 
FOR I31 TO 20 

OUTPUT 709;"AÀS SA" 

Sum=0 

FOR Ja1 TO 5 

ENTER 709;Z 

Sun=Sum+E 

NEXT J 

SavesSum/5 

TempsFNTvsv(Save) 

PRINT TAB(3);I;TAB(15);Temp 
NEXT I 

3EEP 

WAIT 5 

GOTO 590 


PRINT 

OUTPUT 709;"AR AP20 AL22 VRS5" 

POR Is1 TO 3 

OUTPUT 709;"AS SA" 

Sum=0 

FOR J=1 TO S 

ENTER 709; 

Sum=Sum+E 

NEXT J 

IF I=1 THEN Volt=Sum/5 

IF I=2 AND Ido=3 THEN 

PRINT "MAKE SURE VOLTAGE BOX IS SET TO MAIN HEATERS" 
Amp=Sum/ 5 

END IF 

IP I=3 AND Ido=4 THEN 

PRINT "MAKE SURE VOLTAGE BOX IS SET TO AUX HEATERS" 
Amp=Sum/ 5 

END IP 

NEXT I 

Anp=ABS(Amp*1.9182) 

VoltzABS(Volt*25) 

Power=Volt*Amp 

Resistence=Volt/Amp 

PRINT 

BEEP 

PRINT "VOLTAGE(V) CURRENT(A) RESISTENCE( ohms) POWER(W)” 
PRINT 

PRINT USING "1X,5(MDDDD.DD,4X)";Volt,Amp,Resistence,Power 
WAIT S 

GOTO 870 

BEEP 

PRINT 

PRINT "SETUP COMPLETE” 

END 

DEF FNTvsv(V) 

COM /Cc/ C(7) 

T=C(0) 

FOR I=1 TO 7 

ab (ise 


NEXT L 
T=T+8.626897E-2+T*(3.761199B-3-T*5.0689259B-5) 
RETURN T 

PNEND 
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APPENDIX B. PROGRAM DRPJY 


This data acquisition/reduction program is writ nE 
Hewlett-Packard Basic 5.0 language and is listed on the 


following pages. 
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10 
20 
30 
40 
30 
60 
70 
80 
90 
100 
110 
1201 
130 
140 
150 
160 
170 
180 


| FILE NAMB: DRPJY 

| DATE: MAY 15, 1993 

| REVISED VERSION OF DRP72 FOR R-124 (2 INSTRUMENTED TUBES) 
| REVISED BY LT. GEORGE PERRY & LT. JOE YUSICIAN 


COM /Idp/ Idp 

PRINTER IS 1 

CALL Select 

BEEP 

PRINT 

PRINT "DATA COLLECTION COMPLETE" 
END 


SUB Main 

COM /Iprop/ Ift 

COM /Idp/ Idp 

COM /Cc/ C(7),Ical 

COM /Wil/ D2,Di,Do,L,Lu,Kcu 

DIM Em£(20),T(20),D1a(123),D2a(13),D1a(13),Doa(13),La(13),Lua(13),Xcua( 13), 


Et(19),Tn5(4)(15] 


190 DATA 0.10086091,25727.94369, -767345.8295,78025595.81 

200 DATA -9247486589,6.97688E+11,-2.66192E+13,3.94078E+14 

210 READ C(*) 

220 DATA Smooth,High Flux,Turbo-B,High Plux Mod,Turbo-B Mod 

230 READ Tn$(*) 

240 PRINTER IS 701 

250 PRINT 

260 PRINT " Date :",DATES(TIMEDATE) 

270 PRINT 

280 PRINT USING "10X,""NOTE: Program name : DRPJY""" 

290 BEEP 

2T INPUT "ENTER INPUT MODE (0=3054A,1=PILE)",Im 

31 Imo 

3201 INPUT "SELECT FLUID (0sR-114,123R-124)",Ift 

330 Ift=1 

3401 INPUT "1 OR 2 TUBE OPERATION (ENTER 1 OR 2)", Hwmntu 

350 Hwmnntus2 

3601 INPUT "SELECT HEATING MODE (OSELEC; 1=WATER)",Ihm 

370 Ihm=0 

3801 INPUT "ENTER THERMOCOUPLE TYPE (0=NEW, 1=OLD)",Ical 

390 Ical=0 

400 INPUT "GIVE A NAME FOR THE RAW DATA FILE",D2 files 

410 PRINT USING "16X,""New file name: "",15A";D2,fileS 

420 Sizel=20 

430 CREATE BDAT D2 flle$,Sizel 

440 ASSIGN @File2 TO D2_files 

450 ! 

460 1 DUMMY FILE UNTIL Nrun KNOWN 

470 D1 fileSs"DUMMY" 

480 CREATE BDAT D1 file$,Sizel 

490 ASSIGN GFilei TO Di files 

500 OUTPUT QFilel;Dates 

510 BEEP 

520 INPUT “ENTER NR OF DEFECTIVE UPPER TUBE TCS (O#DEFAULT, 2 MAX.)",Idtc 
530 IF Idtc=0 THEN 

540 Ldtc1=0 

550 Ldtc2=0 

560 PRINT USING "16X,""No defective UPPER TUBE TCs exiat""" 

570 END IF 

580 IP Idtc=1 THEN 

590 BEEP 

600 INPUT “ENTER DEFECTIVE UPPER TUBE TC LOCATION (TEMP 1-8)",Ldtcl 
610 PRINT USING "16X,""UPPER TUBE TC ig defective at LOCATION: *",DD";Ldtcl 
620 Ldtc2=0 

630 END IF 

640 IP Idtc=2 THEN 

650 BEEP 

660 INPUT “ENTER DEFECTIVE UPPER TUBE TC LOCATIONS (TEMP 1-8)",Ldtcl,Ldtc2 
670 PRINT USING "16X,""UPPER TUBE TC are defective at LOCATIONS:"",DD,4X,DD";L 
dtc1,Ldtc2 


LS 


680 END IF 

690 IF Idtc»2 THEN 

700 BEEP 

710 PRINTER IS 1 

720 BEEP 

730 PRINT "INVALID ENTRY" 

740 PRINTER IS 701 

750 GOTO 510 

760 END IF 

770 OUTPUT aFilel;Ldtcl,Ldtc2 

780 IF Hwmntu=1 THEN GOTO 1190 

790 INPUT "ENTER NUMBER OF DEFECTIVE LOWER TUBE TCS (O=DEFAULT, 2 MAX)" ,Aidtc 
800 IF Aidtc=0 THEN 

810 Aldtc1=0 

820 Aldtc2=0 

830 PRINT USING "16X,""No defective LOWER TUBE TCs exist""" 

840 END IF 

850 IF Aidtc=1 THEN 

860 BEEP 

870 INPUT "ENTER DEFECTIVE TC LOWER TUBE LOCATIONS (TEMP 9-16)",Aldtcl 
880 PRINT USING "16X,""LOWER TUBE TC is defective at LOCATION:"",DD";Aldtcl 
890 Aldtc2=0 

900 END IF 

910 IF Aidtc=2 THEN 

920 BEEP 

930 INPUT "ENTER DEFECTIVE LOWER TUBE TC LOCATIONS (TEMP 9-16)", ,Aldtc1,Aldtc2 
940 PRINT USING "16X,""LOWER TUBE TC are defective at LOCATIONS:"",DD,4X,DD";A 
ldtci,Aldtc2 

950 END IF 

960 IF Aidtc»2 THEN 

970 BEEP 

980 PRINTER IS 1 

990 BEEP 

1000 PRINT "INVALID ENTRY" 

1010 PRINTER IS 701 

1020 GOTO 790 

1030 END IF 

1040 OUTPUT AFilel;Aldtc1,Aldtc2 

1050 PRINTER IS 1 

1060 BEEP 

1070 PRINT USING "2X,""Select tube number""" 

1080 PRINT USING "6X,""4 Wieland Hard 8 inch""" 

1090 PRINT USING "6X,""5 HIGH FLUX 8 inch""" 

1100! PRINT USING "6X,""6 GEWA-K 40 Fins/in""" 

1110! PRINT USING "6X,""7 GEWA-K 26 Fins/in""" 

1120! PRINT USING "6X,""8 GEWA-T 19 Fins/in""" 

11301 PRINT USING "6X,""9 GEWA-T OR GEWA-TY 26 Fins/in""" 

1140! PRINT USING "6X,""10 THERMOEXCEL-E""" 

11501 PRINT USING "6X,""11 THERMCEXCEL-HE""" 

1160 PRINT USING "6X,""12 TURBO-B""" 

1170! PRINT USING "6X,""13 GEWA-K 19 Fins/in""" 

1180 INPUT Itt 

1190 OUTPUT QFilel;Itt 

1200 PRINTER IS 701 

1210 IF Itts4 THEN PRINT USING "16X,""Tubes: Wieland Hard 8 inch""" 
1220 IF Itts12 THEN PRINT USING "16X,""Tubes:  Turbo-B""" 

1230 IF Itt=5 THEN PRINT USING "16X,""Tubes: High Flux 8 inch""" 
1240! IF Itt<10 THEN PRINT USING "16X,""Tube Number: e DATE? 
12501 IF Itt»9 THEN PRINT USING "16X,""Tube Number: nn, DD" TEO 
1260! INPUT “ENTER OUTPUT VERSION (0=LONG, 12SHORT, 2=NONE)",Iov 
1270 Iov=0 

1280! INPUT “SELECT (0=LIQ,1=VAP,2=(LIQ+VAP)/2)",Ilqv 

1290 Ilqv=2 

1300! 

1310! DIMENSIONS FOR TESTED TUBES 

13201 

1330! Di=Diameter at thermocouple positions 

1340 DATA .0111125,.0111125,.0111125,.0129540,.012446,.0129540,.0100965 
1350 .DATÀ .01009655,.01157,.01197,.01157, 017157, 01157, 50100965 
1360 READ Dla(*) 

1370 Di=Dla(Itt) 
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13801 

1390! D2=Diameter of test section to the base of fins 

II0SNUDATAT-015875. 0158757 0158757 015824 ..015875,.015824, .01270 
SUE 70127,.0138,.0138,.0138,.0138, .0138,.0127 

1420 READ D2a(*) 

1430 D2=D2a( Itt) 

1440! 

1450! Di=Inside diameter of unenhanced ends 

INR DA EAN OZ 7.0127, .0127,.0132, .0127,.0132,,0111125,:.0111125 
iq QI DATAW 0118, 0118, .0118, .0118, -0118, 0111125 

1480 READ Dia(*) 

1490 Di=Dia(Itt) 

15001 

1510! Do=Outside diameter of unenhanced ends 

RIR UATA 2015875,.015875,.015875,.01582477015875,.015924,.01270,.01270 
Mn ODATA :01331,.01331,.01331,.01331,.0158,.0127 

1540 READ Doa(*) 

1550 Do=Doa(Itt) 

15601 

15701 L=Length of enhanced surface 

IESDESDATAS: 1016, .1015,.10165,.1015,.2032,.2032,.2032,.2032,.2032, .20312,.2032, .20 
H2 2032... 2032 

1590 READ La(*) 

1600 L=La(Itt) 

16101 

1620! Lu=Length of unenhanced surface at the ends 

BS DATAS 0254,.0254,.0254,.0254,.0762,.0762,.0762,.0762,.0762,.0762,.0762,.07 
6027. 0762,.0762 

1640 READ Lua(*) 

1650 Lu=Lua(Itt) 

16601 

16701 Kcu=Thermal Conductivity of tube 

1680 DATA 398,344,344,45,344,45,344,344,398,398,398,398,398,344 
1690 READ Kcua(*) 

1700 Kcu=Kcua(Itt) 


17101 
1720 Xn=.8 
1730 Fr=.3 


1740 C€=3.7037E+10 

1750 A=PI*(Do*2-Di~2)/4 

1760 P=PI*Do 

1770 Nsub=0 

1780 Ntc=20 

1790 J=1 

1800 Repeat: | 

1810 Dtld=2.22 

1820 Ido=2 

1830 ON KEY 1,15 RECOVER 1800 

1840 PRINTER IS 1 

1850 PRINT USING "2X,""SELECT OPTION""" 

1860 PRINT USING "6X,""O-TAKE DATA""" 

1870 PRINT USING "6X,""1=SET UPPER TUBE HEAT FLUX (MAIN)""" 
1880 PRINT USING "6X,""2=SET Tsat""" 

1890 PRINT USING "6X,""3=SET LOWER TUBE HEAT FLUX (AUX)""" 
1900 BEEP 

1910 INPUT Ido 

1920 IF Ido>3 THEN Ido=3 

1930 IF Ido=0 THEN 3680 

19401 

19501 LOOP TO SET UPPER TUBE HEAT FLUX 

1960 IF Idos1 THEN 

1970 OUTPUT 709;"AR AF20 AL21 VR5" 

1980 PRINT "SET VOLT BOX TO MAIN" 

1990 DEER 

2000 INPUT "ENTER DESIRED Qdp",Dqdp 

2010 PRINT USING "4X,""DESIRED Qdp ACTUAL Qdp""" 


ld 


2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
23701 
23801 
2390 
2400 
2410 
24201 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 


Err=1000 

FOR Is1 TO 2 

OUTPUT 709;"AS SA" 
Sum=0 

FOR Ji=1 TO 5 

ENTER 709;E 

Sum=Sum+E 

NEXT J1 

IF I=1 THEN Volt=Sum/5 
IF I=2 THEN Amp=Sum/5 
NEXT I 
Amp-ABS(Amp*1.9182) 
Volt=ABS(Volt*25) 
AqdpsVolt^Amp/(PI*D2*L) 
IF ABS(Aqdp-Dqdp)>Err THEN 
IF Aqdp»Dqdp THEN 

BEEP 4000,.2 

BEEP 4000,.2 

BEEP 4000,.2 

ELSE 

BEEP 250,.2 

BEEP 250,.2 

BEEP 250,.2 

END IF 

PRINT USING "4X,MZ.3DE, 2X,MZ.3DE";Dqdp,Aqdp 
WAIT 2 

GOTO 2030 

ELSE 

BEEP 

PRINT USING "4X,MZ.3DE,2X,MZ.3DE";Dqdp, Aqdp 
Err=500 

WAIT 2 

GOTO 2030 

END IF 

END IF 


LOOP TO SET Tsat 

IP Ido=2 THEN 

IF Ikdt=1 THEN 2440 
Dtld=2.22 

PRINT USING "4X,"" DTsat ATsat Rate Tv Rate or 
Ikdt=1 

Old120 

01d2=0 

Nn=1 

Nrs=Nn MOD 15 

NnsNn*1 

IP Nrs=1 THEN 

PRINT USING "4X,"" Tsat Tdi Tld2 Tv Tsunmp""" 
END IF 

OUTPUT 709;"AR AF16 AL19 VR5" 
FOR I=1 TO 6 

IF I>4 THEN 2700 

Sum=0 

OUTPUT 709;"AS SA" 

FOR Jiz1 TO 20 

ENTER 709;Eliq 
Sum=Sum+Eliq 

NEXT J1 

Eliq=Sum/20 
TldsFNTvsv(Eliq) 

IP I=1 THEN Tldl-Tld 
IP I=2 THEN Tld2=Tld 
IF I*3 THEN Tv4Tld 

IP 1-4 THEN Tsump=T1d 
IP I=5 THEN Tinlet=Tld 
IF I=6 THEN Tout=Tld 


TRO 


2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
31801 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 


NEXT I 

Atld=(Tldl+Tld2)".5 

IF ABS(Atld-Dtld)>.2 THEN 

IF Atld>Dtld THEN 

BEEP 4000,.2 

BEEP 4000,.2 

BEEP 4000,.2 

ELSE 

BEEP 250,.2 

BEEP 250,.2 

BEEP 250,.2 

END IF 

Errl=Atld-Oldl 

Oldl-Atld 

Err2=Tv-Old2 

Old2=Tv 

IF T1d1>100. THEN 2880 

PRINT USING "4X,5(MDD.DD,2X)";Dtld,Tld1,Tld2,Tv,Tsump 
WAIT 2 

GOTO 2470 

ELSE 

IF ABS(Atld-Dtld)».1 THEN 

IF Atld»Dtld THEN 

BEEP 3000,.2 

BEEP 3000,.2 

ELSE 

BEEP 800,.2 

BEEP 800,.2 

END IF 

Errl=Atld-Oldl 

Oldi=Atld 

Err2=Tv-Old2 

Old2=Tv 

PRINT USING "4X,5(MDDD.DD,2X)";Dtld,Tld1,T1Q2,Tv,Tsump 
WAIT 2 

GOTO 2470 

ELSE 

BEEP 

Errl=Atld-Oldl 

Oldi=Atld 

Err2=Tv-Old2 

Old2=Tv 

PRINT USING "4X,5(MDDD.DD,2X)";Dtld,Tld1,Tld2,Tv,Tsump 
WAIT 2 

GOTO 2470 

END IF 

END IF 

END IF 

l 

LOOP TO SET LOWER TUBE HEAT FLUX 
IF Ido=3 THEN 

PRINT " SET VOLT BOX TO AUX" 
OUTPUT 709;"AR AF20 AL22 VRS" 
BEEP 

INPUT "ENTER DESIRED AuQdp",Duxqdp 
PRINT USING "2X,""DESIRED AuxQdp ACTUAL AuxQdp""" 
Errz1000 

FOR I=1 TO 3 

OUTPUT 709;"AS SA" 

Sum=0 

FOR Ji=1 TO 5 

ENTER 709;E 

Sum=Sum+E 

NEXT Ji 

IF I=1 TEEN Volt=Sum/5 

IF I=3 THEN Amp=Sum/5 

NEXT I 


ss 


3360 Amp=ABS(Amp*1.9182) 

3370 Volt=ABS(Volt*25) 

3380 Auxgqdp=Volt*Amp/(PI*D2*L) 

3390 IF ABS(Auxqdp-Duxgqdp)>Err THEN 
3400 IF Auxqdp»Duxqdp THEN 

3410 BEEP 4000,.2 

3420 BEEP 4000,.2 

3430 BEEP 4000,.2 

3440 ELSE 

3450 BEEP 250,.2 

3460 BEEP 250,.2 

3470 BEEP 250,.2 

3480 END IF 

3490 PRINT USING "4X,MZ.3DE,2X,MZ.JDE";Duxqdp,Auxqdp 
3500 WAIT 2 

3510 GOTO 3260 

3520 ELSE 

3530 BEEP 

3540 PRINT USING "4X,M2.3DE,2X,MZ.3DE";Duxqdp,Auxqdp 
3550 Err=500 

3560 WAIT 2 

3570 GOTO 3260 

3580 END IF 

3590 GOTO 1850 

3600 END IF 

3610! ERROR TRAP FOR Ido OUT OF BOUNDS 
3620 IF Ido>3 THEN 

3630 BEEP 

3640 GOTO 1850 

3650 END IF 

36601 

3670! TAKE DATA LOOP 

3680 IF Ikol=1 THEN 3730 

3690 BEEP 

3700 Ikol=1 

3710 INPUT “ENTER BULK OIL %",Bop 
3720 Ikol=1 

3730 OUTPUT 709;"AR AFOO AL19 VR5" 
3740 Ntc=20 

3750 FOR I=1 TO Ntc 

3760 OUTPUT 709;"AS SA" 

3770 Sum=0 

3780 FOR Ji=1 TO 20 

3790 ENTER 709;E 

3800 Sum=Sum+E 

3810 IF I=(17-Nsub) OR I=(18-Nsub) THEN Et(J1-1)=E 
3820 NEXT Ji 

3830 Kdl=0 

3840 IF I=(17-Nsub) OR I=(18-Nsub) THEN 
3850 Eave=Sum/20 

3860 Sum=0. 

3870 FOR Jk=0 TO 19 

3880 IF ABS(Et(Jk)-Eave)<5.0E-6 THEN 
3890 Sum=Sum+Et(Jk) 

3900 ELSE 

3910 KdlsKdl*1 

3920 END IF 

3930 NEXT Jk 

3940 IF I=(17-Nsub) OR I=(18-Nsub) THEN PRINT USING "4X,""Kdl = "",DD";Kdl 
3950 IF Kdl>20 THEN 

3960  BEEP 

3970  BEEP 

3980 PRINT USING "4X,""Too much scattering in data - repeat data set""" 
3990 GOTO 1840 

4000 END IF 

4010 END IP 

4020 Emf(I)=Sum/(20-Kdl) 
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4030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290! 
43001 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4490 


NEXT I 

Coun=0. 

OUTPUT 709;"AR AF20 AL22 VR5" 
FOR I=1 TO 3 

OUTPUT 709;"AS SA" 

Sum=0 

FOR Ji=1 To 5 

ENTER 709;E 

Sum=Sum+E 

NEXT J1 

IF Coun=0. THEN 

IF I=1 THEN Vr=Sum/5 

IF I=2 THEN Ir=Sum/5 

ELSE 

IF I=1 THEN Avr=Sum/5 

IF I=3 THEN Air=Sum/5 

END IF 

NEXT I 

IF Coun=0 THEN 

PRINT “SHIFT VOLT BOX TO AUX" 
END IF 

IF Coun=1 THEN GOTO 4290 
INPUT "TAKE AUX READINGS(1=YES)?",Ccon 
IF Ccon>1 THEN GOTO 4250 
Coun=Coun+Ccon 

GOTO 4060 


CONVERT emf'S TO TEMP,VOLT,CURRENT 
Twa=0 

Atwa=0 

FOR I=1 TO Ntc 

IF Idtc>0 THEN 

IF I=Ldtcl OR I=Ldtc2 THEN 
T(I)s-99.99 

GOTO 41470 

END IF 

END IF 

IF Aidtc>0 THEN 

IF I=Aldtci OR I=Aldtc2 THEN 
T(I)=-99.99 

GOTO 4470 

END IF 

END IF 

T(I)=FNTvsv(Emf(I)) 

NEXT I 

Twa=Twa/(4-Idtc) 

FOR I=1 TO 8 

IF I=Ldtcl OR I=Ldtc2 THEN 
Twa=Twa 

ELSE 

TwazTwa+T(I) 

END IF 

NEXT I 

Tw=Twa/(8-Idtc) 

FOR I=9 TO 16 

IF I=Aldtcl OR I=Aldtc2 THEN 
Atwa=Atwa 

ELSE 

Atwa=AtwatT(I) 

END IF 

NEXT I 

Atw-Atwa/(8-Aidtc) 
Tld=T(17-Nsub) 
Tld2=T(18-Nsub) 
Tida=(Tld+T1ld2)*.5 
TvzT(19-Naub) 
TaumpsT(20-Nsub) 


141 


4700 
4710 
4720 
4730 
4740 
4750 
4760 
47701 
47801 
4790 
4800 
4810 
48201 
4830! 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980! 
49901 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
51101 
5120: 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 
5210 
52201 
52301 
52401 
5250 
5260 
5270 
5280 
5290 
5300 
5310 
5320 
5330 
5340 
5350 
53601 
53701 
53801 
5390 
5400 
5410 


AmpsABS(Ir*1.9182) 
Volt-ABS(Vr)*25 
Q-Volt*Amp 
Auamp=ABS(Air*1.9182) 
Auvolt-ABS(Avr)*25 
Auq=Auvolt*Auamp 
KcuszKcua(Itt) 


FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
Tw=Tw-Q*LOG(D2/D1)/(2*PI*Kcu*L) 
Atw=Atw-Auq*LOG(D2/D1)/(2*PI*Kcu*L) 

Tsat-(Tlda*Tv)*.5 


COMPUTE THETHAB 

Thetab-Tw-Tsat 

Athetab=Atw-Tsat 

IF Thetab<O THEN 

BEEP 

INPUT "UPPER TWALL<TSAT (O=CONTINUE, 1-END)",Iev 
IF Iev=0 THEN GOTO 3150 

IF Iev=1 THEN 1850 

END IF 

IF Athetab<0O THEN 

BEEP 

INPUT "LOWER TWALL<TSAT (O=CONTINUE, 1=END)",Aiev 
IF Alev=0 THEN GOTO 3150 

IF Alev=1 THEN 1850 

END IF 


COMPUTE VARIOUS PROPERTIES 
Tfilm=(Tw+Tsat)*.5 
Rho-FNRho(Tfilm) 
Mu=FNMu(Tfilm) 
K=FNK(Tfilm) 
Cp=FNCp(Tfilm) 
Beta=FNBeta(Tfilm) 
Hfg-FNHfg(Tsat) 
Ni-Mu/Rho 
Alpha=K/(Rho*Cp) 
Pr=Ni/Alpha 
Psat=FNPsat(Tsat) 


COMPUTE VARIOUS PROPERTIES FOR AUX TUBE 
Atfilm=(Atw+Tsat)*.5 

Arho=FNRho(Atfilm) 

Amu=FNMu(Atfilm) 

Ak=FNK(Atfilm) 

Acp=FNCp(Atfilm) 

Abeta=FNBeta(Atfilm) 

Ani=Amu/Arho 

Aalpha=Ak/(Arho*Acp) 

Apr=Ani/Aalpha 


COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
FOR UNENHANCED END(S) 

Hbar=190 

Fe=(Hbar*P/(Kcu*A))”.5*Lu 

TanhzFNTanh(Fe) 

Theta-Thetab*Tanh/Fe 
Xx=(9.81*Beta*Thetab*Do”3*Tanh/(Fe*Ni*Alpha))”.166667 
Yy-2(1*(.559/Pr)^(9/16))^(8/27) 
Hbarc-K/Do*(.6*.387*Xx/Yy)^2 

IF ABS((Hbar-Hbarc)/Hbarc)».001 THEN 
Hbars(Hbar*Hbarc)*.5 

GOTO 5260 

END IF 


COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
FOR UNENHANCED END(S) FOR AUX TUBE 

Ahbar=190 

Fez(Ahbar*P/(Kcu*A))^.5*Lu 

Atanh=FNTanh(Fe) 
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5420 Atheta=Athetab*Atanh/Fe 

5430 Axx=(9.81*Abeta*Athetab*Do~3*Atanh/ (Fe*Ani*Aalpha) )~*.166667 
5440 Ayy=(1+(.559/Apr)*(9/16) )*(8/27) 

5450 Ahbarc=Ak/Do*(.6+.387*Axx/Ayy)72 

5460 IF ABS( (Ahbar-Ahbarc)/Ahbarc)>.001 THEN 

5470 Ahbar=(Ahbar+rAhbarc)*.5 

5480 GOTO 5400 

5490 END IF 

5500! 

5510! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) 

5520 Q1=(Hbar*P*Kcu*A)*.5*Thetab*Tanh 

5530 Qc=Q-2*Q1 

5540 As=PI*D2*L 

5550 ! 

5560 !COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) OF AUX TUBE 
5570 Aql=(Anbar*P*Kcu*A)”.5*Athetab*Atanh 

5580  Aqc-Auq-2*Aql 

55901 

56001 COMPUTE ACTUAL HEAT FLUX AND BOILING COEFFICIENT 

5610 Qdp=Qc/As 

5620 Htube=Qdp/Thetab 

5630 Csf-(Cp*Thetab/Hfg)/(Qdp/(Mu*Hfg)*(.014/(9.81*Rho)^.5)^(1/3. ) *Pr^1.7) 
56401 

56501 COMPUTE ACTUAL HEAT FLUX AND BOILING COEFFICIENT FOR AUX TUBE 
5660  AuqdpsAqc/As 

5670 Ahtube=Auqdp/Athetab 

5680 Acsf-(Acp*Athetab/Hfg)/(Auqdp/(Amu*Hfg)*(.014/(9.81*Arho)^.5)^(1/3.) *Apr^1 
7) 

56901 

57001 COMPUTE CHURCHILL-CHU CORRELATION FOR UPPER TUBE 

5710 Cc=Bbar*Thetab 

57201 

57301 RECORD TIME OF DATA TAKING 

5740 OUTPUT 709;"TD" 

5750 ENTER 709;ToldS 

57601 

57701 OUTPUT DATA TO PRINTER 

5780 PRINTER IS 701 

5790 PRINT 


5800 PRINT USING "10X,""Data Set Number - "",DDD,2X,""Bulk Oil $ - "",DD.D";J,B 

op 

5810 PRINT " TIME:",TIMES(TIMEDATE) 

5820 PRINT 

5830 PRINT USING "10X,""TC No: 1 2 B 4 5 6 7 
guo 


5840 PRINT USING "10X,""Temp :”"”,8(1X,MDD.DD)";T(1),T(2),T(3),T(4),T(5),T(6),T( 

7),T(8) 

5850 PRINT USING "10X,""TC No: 9 10 11 12 13 14 15 
16""" 

5660 PRINT USING "10X,""Temp :"",8(1IX,MDD.DD)";T(9),T(10), T( 11); TC 12) , T( 13) , T(1 

4),T(15),T(16) 

5870 PRINT 

5880 PRINT USING "10X,"" Twa  ATwa  Tliqd  Tliqd2 Tvapr Psat  Tsat “su 

m mem 

5890 PRINT USING "10X, 3(MDD.DD, 1X), 1X,MDD.DD, 1X, 3( 1X,MDD.DD),2X,MDD.D";Tw,Atw,T 

ld,Tld2,Tv,Psat,Tsat,Tsump 


5900 PRINT USING "10x,"" Vr Ir ES 
5910 PRINT USING "10X,2(MDDD.DD,3X),1X,MZ.3DE";Vr,Ir,CC 
5920 PRINT USING "10X,"" Thetab Htube Qdp Athetab  Ahtube AuQd 


5930 PRINT USING "10X,MDD.3D, 1X,MZ.3DE, 1X,MZ. 3DE, 1X, MDD. 3D, 1X,MZ. 3DE, 1X,MZ. 3DE”" 
; Thetab, Htube, Qdp, Athetab, Ahtube, Auqdp 

5940 PRINT 

5950 BEEP 

5960 INPUT "OK TO STORE THIS DATA SET (1=Y,0=N)?",Ok 
5970 IF Ok=1 THEN J=J+1 

5980 IF Ok=1 AND Im=0 THEN 

5990 OUTPUT GFilel;Bop,Told$,Emf(^*),Vr,Ir,Avr,Air 
6000 END IF 

6010 IF Iuf=1 THEN OUTPUT QUfile;Vw,Uo 

6020 IF Ire=1 THEN OUTPUT @Refile;Fms, Rew 

6030 DEER 
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6040 INPUT "WILL THERE BE ANOTHER RUN (1-Y,0-N)?",Go on 
6050 Nrun=J 

6060 IF Go. on30 THEN 6080 

6070 IF Go_on<>0 THEN Repeat 

6080 BEEP 

6090 PRINT 

6100 PRINT USING "10X,""NOTE: "",Z22,"" data runs were stored in file "",15A";J- 
1,D2 filos 

6110 ASSIGN @Filel TO * 

6120 OUTPUT @File2;Nrun-1 

6130 ASSIGN QGFilel TO D1 fileS 

6140 ENTER @Filel;DateS,Ldtcl,Ldtc2,Itt 
6150 ENTER @Filel;Aldtcl,Aldtc2 

6160 OUTPUT @File2;DateS,Ldtcl,Ldtc2,Itt 
6170 OUTPUT QFi1le2;Aldtcl,Aldtc2 

6180 FOR I=1 TO Nrun-1 

6190 ENTER SFilel;Bop,ToldS,Emf(*),Vr,Ir 
6200 ENTER GFilel;Avr,Air 

6210 OUTPUT @File2;Bop,Told$,Emf(*),Vr,Ir 
6220 OUTPUT @File2;Avr,Air 

6230 NEXT I 

6240 ASSIGN @Filel TO * 

6250 PRINTER IS 1 

6260 PURGE "DUMMY" 

6270 PRINT 

6280 PRINT "DUMMY FILE PURGED" 

6290 SUBEND 

63001 

63101 CURVE FITS OF PROPERTY FUNCTIONS 
63201 

6330 DEF FNKcu(T) 

63401 OFHC COPPER 250 TO 300 K 

6350 Tk=T+273.15 [C TO K 

6360 K=434-.112*Tk 

6370 RETURN K 

6380 FNEND 

6390! 

6400 DEF FNMu(T) 

6410 COM /Iprop/ Ift 

6420 IF Ift=1 THEN 

6430 1223 TO 373 K 

64401 170 TO 360 K CURVE FIT OF VISCOUSITY 
6450 Tk=T+273.15 E TO sm 

6460 MusEXP(-4.4636*(1011.47/Tk)) *1.0E-3 
6470 END IF 

6480 IF Ift=1 THEN 

6490! 223 TO 373 K 

6500 Muz23.4593-4.26E-2*T*3.9485E-4*T^2-4.193E-6*T^3*2.0709E-8*T^4 
6510 Mu=Mu*1.0E-4 

6520 END IF 

6530 RETURN Mu 

6540  FNEND 

65501 

6560 DEF FNCp(T) 

6570 COM /Iprop/ Ift 

6580 IF Ift-O THEN 

65901 180 TO 400 K CURVE FIT OF Cp 

6600 TksT*273.15 IC TO K 

6610 Cpx=.40188+1.65007E-3*Tk+1.51494E-6*Tk"2-6.67853E-10*Tk” 3 
6620 Cp=Cp*1000 

6630 END IF 

6640 IF Ift=1 THEN 

6650! 223 TO 373 R CURVE FIT Cp(R124) 

6660 Cps1.0542*2.1405E-3*T*1.0709E-5*T^2-6.4721E-8*T^3-1.4324E-9*T^4* 4. 136E-11* 
T^ 

6670 Cp=Cp*1000 

6680 END IF 

6690 RETURN Cp 

6700 FNEND 


144 


67101 

6720 DEF FNRho(T) 

6730 COM /Iprop/ Ift 

6740 IF Ift-O THEN 

5750 Tk=T+273.15 {Cr TOLER 

6760 X=1-(1.8*Tk/753.95) IK TO R 

6770 Ro=36.32+61.146414*X^(1/3)+16.418015*X+17.476838*X^.5+1.119828*X^2 
6780 Ro=Ro/.062428 

6790 END IF 

6800 IF Ift=1 THEN 

6810! 223 TO 373 K (R124) 

6820 Ro=1434.8-2.8619*T-6.7267E-3*T^2-7.2852E-5*T^3 
6830 END IF 

6840 RETURN Ro 

6850 FNEND 

6860! 

6870 DEF FNK(T) 

6880 COM /Iprop/ Ift 

6890 IF Ift=0 THEN 

6900! T<360 K WITH T IN C 

6910 K=.071-.000261°T 

6920 END IF 

6930! RETURN K 

6940 IF Ift=1 THEN 

6950! 223 TO 373 K (R-124) 

6960 R=7.5191E-2-3.5436E-4*T-1.9545E-7*T"2+3.1835E-9*T"J3 
6970 END IF 

6980 RETURN K 

6990 FNEND 

70001 

7010 DEF FNTanh(X) 

7020 P=EXP(X) 

7030 Q=1/P 

7040 Tanh=(P-Q)/(P+Q) 

7050 RETURN Tanh 

7060 FNEND 

7070! 

7080 DEF FNTvsv(V) 

7090 COM /Cc/ C(7),Ical 

7100 T=C(0) 

7110 FOR I=1 TO 7 

DIU TeT«C(T)*V^I 

7130 NEXT I 

7140 IF Ical=1 THEN 

7150 T=T-6.7422934E-2+T*(9.0277043E-3-T*(-9.3259917E-5)) 
7160 ELSE 

7170 T=T+8.626897E-2+T*(3.76199E-3-T*5.0689259E-5) 
7180 END IF 

7190 RETURN T 

7200 FNEND 

me LO! 

7220 DEF FNBeta(T) 

7230 Rop=FNRho(T+.1) 

7240 Rom=FNRho(T-.1) 

7250 Beta=-2/(Rop+Rom)*(Rop-Rom)/.2 

7260 RETURN Beta 

7270 FNEND 

72801 

7290 DEF FNHfg(T) 

7300 COM /Iprop/ Ift 

7310 IF Ift=0 THEN 

7320 Hfg=1.3741344E+5-T*(3.3094361E+2+T*1.2165143) 
7330 END IF 

7340 IF Ift=1 THEN 

7350! 223 TO 373 K CURVE FIT (R124) 

7360 Hfg=159.93-.4609*T-1.458E-3*T"2-1.3715E-5*T"3 
7370 Hfg=Hfg*1000 
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7380 
7390 
7400 
7410! 
7420 
7430 
7440 
7450! 
7460 
7470 
7480 
7490 
7300 
7510! 
7520 
7530 
7540 
7550 
7360 
7370 
7580 
7590 
7600 
7610 
7620 
7630 
7640 
7650 
7660 
7670 
7680 
7690 
7700 
7710 
7720 
7730 
7740 
7750 
7760 
7770 
7780 
7790 
7800 
7810 
7820 


END IF 
RETURN Hfg 
FNEND 


DEF SNPsat(Tc) 

COM /Iprop/ Ift 

IF Ift=0 THEN 

Y TO BO deg = CURVE ZIT CF 2335 

Tis].3*Tct32 
Pa-5.945525*75*(.15352082-TT*(1.194096352-2* 0549. 515057 28 $500) 
PgsPa-14.7 

END IF 

IF Ift=1 THEN 

223 TO 373 X CURVE SIT (R-124) 
Pa=(162.68+5.946*Tc+9.2081E-2*Tc"2+6.9359E-4*Tc0c"3)*.14504 
2g=Pa-14.7 

END IF 

IF Pg»0 THEN | +aPSIG,-=in dg 

PsatsPg 

ELSE 

Psat=Pg*29.92/14.7 

SND IF 

RETURN Psat 

FNEND 


46. ¿ 


SUB Purg 
3EEP 

INPUT "ENTER FILE NAME TO BE DELETED", Files 
PURGE Files 

GOTO 7650 

SUBEND 


SUB Select 

COM /Idp/ Idp 

BEEP 

PRINTER IS 1 

PRINT USING "2X,""Select option:""" 
PRINT USING "6X,"" O TAKE DATA""" 
PRINT USING "SX,"" 1 PURGE ?7ILES""" 
INPUT Idp 

IF Idp=0 THEN CALL Main 

IF Idp=1 THEN CALL 2urg 

SUBEND 
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APPENDIX C. 


REPRESENTATIVE DATA SET 


Jate 7 Seg 1395 
NOTE: Program name : ORPJY 
New file name: 009075C2 
UPPER TUBE TC ara defactive at LOCATIONS: 7 2 
SD ELL T at LOCATION: 9 
Tubes: Wisland Hard 3 incn 
Data Set Number = | sulk Orie 2.4 
TIME: 15515: 58 
TC No: | E 3 4 3 9 7 g 
T amp 3.8S 303 Jus 3.56 Ses 3.39 -39.39 -39.39 
TC No: 9 18 11 Ee 13 14 IS 18 
Temp : -99.39 29.29 18.36 292.46 19.604 20.49 13.53 17.30 
Twa ATua Tliad Tliad2 Tvapr Psat Tsat Tsump 
3:29 19.36 20790 eee sz TIS? 2.18 -10.3 
Ve Ir Dc 
.34 Se 2. 388E+02 
Thetan Htube Qdo Athetab Ahtube out do 
1.525 3.48088E+02 S.342£+02 17.194 1.832€-92 3.150E-23 
Qata Set Number = 2 Bulk 011 0.0 
TIME: {3:14:38 
TC No: t Z 3 4 S 6 7 3 
Tamp 4.37 4.31 de de 2236 4.98 S19 239 390239 39 
TC No: q 19 RI I l3 } 4 1S 16 
Tepo LO pal 29:57 20. 07 29.32 19.383 17.58 
Twa ATwa Tliad Tliad2 Tvapr Psat Tsat Tsump 
4,94 19.55 EE Zou 2.18 19.30 2.28 -10.8 
Ur ER Cc 
pus seo 4.530E-02 
Thetan Htube Qdo Athetab | Ahtube AuQdgo 
2.538 3.363E€+02 .1.340E-03 17.397 1.793E+02 3.120€+03 
Data Set Number = 3 Bulk O11 % = 9.2 
INES 13:16:14 
TO No: 1 2 5 4 S S 7 3 
Temp 7.38 god T 6.32 ug ds asas 13.39 
TO Nos 3 18 11 lez 13 14 1S 16 
Temp 2423.39 790.08 018775 0020.32 2203541. -20.18- 18.35 16.96 
Tua ATua Ie ad Tliad2 Tvaor Psat Tsat Tsump 
S. T2 0% Za 2-025. 509573 2.08 -10.3 
Ur Dr Ee 
1.49 Ge 1. I2SE+O3 
Thetab Htube Jap Athetab Ahtude Ou Op 
3.104 3.41!9E-«02 1.740E+93 17.957 1.360€+02  3.!72E*903 
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Qata 3e 
TIME: 


TC No: 
Tamo 
nC Mo: 
Tema 


Tua 
13.11 
yr 
Eë 
Thetao 
SPIRA 


Qata Se 
TIME: 


TC No: 
Temp 

TE No: 
Tamo : 


Tua 

14.95 
Ur 
NS 
Thetao 
l! .362 


Data Se 
TIME: 


TC No: 
Temo 
TC Na: 
Tema 


Tua 

| 07548 
Ur 
cud 
Thetab 
15.2324 


t Number = 
JS. Saa 
i E 
19.10 19.19 
2 19 
=39 ag E da 
ATwa ened 
19.18 Ss 
Ir 
233 
Htube 
3.319€+92 
t Number = 
13:18:54 
| 2 
Acao SIR 
g 10 
=99.99 20.89 
ATua Tliad 
ES 2.91 
Ir 
.39 
Htube 
3: :550E r02 
t Number = 
13:29:03 
| Z 
PE Nr AP L^ 
g 19 
IZ 
ATua lad 
DS te 2.94 
Ir 
as 
Htube 


3.37939E+02 TS OATES OS 


4 


- 
Ze 


S 


4 


3ulk O11 % = 3.3 
5 d 5 5 i 3 
POSE 3.64 19.91 10.zð -329.29 -39.39 
L1 12 13 14 1S 15 
[8.77 29.12. 19:31 20.29 A co 
TlLrad2 Tvapr Psat Tsat Tsumg 
2.18 1.92 19.36 1.99 -19.7 
Ec 
1.895E+03 
Qdo Athetao Ahtube AâuQdo 
.395€+03 17.189 1.334E-02 3.1928+03 
Bulk Qil X = Q.2 
3 4 S 5 7 3 
14.12 13.47 13.7S 14.50 -39.99 -99.39 
11 E ES 14 15 18 
19.23 d 20.92 223.32 918-37 las 
Tliad2 Tvapr Psat T3at Tsuma 
2.28 2.04 10.76 2.09 -10.8 
cc 
3.3468-93 
Qdo Athetab Ahtube AuQda 
.342£r03 17.429 1.790E+92 3.102£-03 
Bulk Qil X > Q.2 
3 4 3 S 7 3 
17.30 18.31 16.54 17.53 -99.9339 -99.99 
11 D iS 14 IS 18 
18.88 20.32 13.383 20.48 18.47 17.38 
Tl:gd2 Tvapr Psat Tsat Tsumo 
A oie 2.39 19.30 2 0 32 7 
Gc 
4.984E+93 
Qdn Athetab Ahtube AuQdg 


17,139 1.334E+02 3.1435E+03 
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Data Set 
TIME: 


=o Ne: 
Tamo 
TERNO: 


Number = 
Ier E 5 
| 2 
[9.49 13.10 
g 19 


Mano 99. 39. 393-33 


Twa 
gg 
Ur 
2.80 
Thetao 
16.783 


Data Set 
TIME: 


TC No: 
Temp 
TC No: 


ATwa rada 


13:23 as 
Ir 
.50 
Htube 


5.312 2 


Number s 
iss] 


I 2 
E OO peca 
g 19 


Temo : -99.99 19.75 


Tua 
21.34 
Ur 
agr 
Thetau 
19.389 


ATwa Tliqd 

18.37 1.31 
Ir 
.85 


Data Set Number = 


TIME: 


TE No: 
Temp 
TC No: 


Temp : -99.99 20. 


Twa 
26.93 
Ur 
Sus 
Thetab 
23.797 


ee 

I 2 
ce Jl. (ono) 
3 19 
33 


ATwa Tliad 
139.59 Bold 
ic 
. 79 
Htube 


7 Bulk Jl > = 2.92 


6 


8 


3 


3 4 5 5 7 3 
13.35 18.05 13.12 13.37 -29.39 -39.39 
[1 e I3 14 15 [8 
IS 9. 34 909.350 —30.11 13.08 17.07 

Tligd2 Tvapr Psat Tsat Tsuma 
gl 1.38 19.58 2.90 -19.7 
Sc 

4.707£+03 
Qdo Athetab Ahtube Audo 


:d5lE*J3 17.02% —1.848b6f02 3.14595E+03 


Bulk Oil 2 = 29.9 
3 4 S 5 7 3 
Zn. 96 Meu see Da 22 a +99. 39 «99 39 
E 12 13 14 1S 16 


13.33 6 13:73 19.35 20.98 18.904 (18.35 


Tliad2 Tvapr Psat Tsat Tsump 


219 1-89 10.63 1,98 -10.7 
ES 

S.3S08+03 

Qdo Athetab Ahtube AuQcp 


2495&7033 250329 VI. GBMESDZ 5,1526403 


Bulk Oil à = 98.23 
5 4 5 5 7 g 
29 UNE S 1208243290079 5.001 3391339 89 39.39 
11 E ES 14 IS 16 


19. 2520. colo Dadas. ess 17.71 


Tliad2 Tvapr Psat Tsat Tsuma 
2.47 2093 NEEN D CT Sam 
Cc 

7.3639£+03 
Qdg Athetab  Ahtube ñuQdp 


3.391£+92 3.498586E«*03 17.328 1.78S5E*02  Á3.094E*03 
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Dara Set Numoar = 


TIME: 


TOIT NO: 
Tamo 
TC No: 


bo er 028 


| E 
lG 11.52 
3 19 


Tamo : -99.39 17.56 


Tua 
BEL 
Ur 
Dong 
Thetab 
30931 


Data Set Number = 


TIME: 


TC No: 
Temp 
TC No: 
Temp 


Tua 
12.14 
Ur 
4.13 
Thetab 
19.927 


Qata 3e 
TIME: 


TC No: 
Temg : 
TC No: 
Temg 


Twa 
14.37 
Ur 
5.38 
Thetab 
2859 


ATwa Titad 
16.34 EE 
Ir 
s 
Htube 


1.310605 


lo 2g 
| 2 
LECHE Ss 
3 19 
ISO 
ATwa Mia 
ara 2.23 
Ir 
. 356 
Htube 
532 218 053 


t Numoer = 12 


13:30:20 

l E 
13.30 14.46 
g 19 
S miGec! 
ATua Tiliga 
i S 2.1 

lr 

1.25 

Htube 


12 


E 


11 


E 


Julk Orla" 0.23 


S 4 E 3 » 2 
no IP cU SA 355535 555 3 5) 
11 Pz 13 14 15 18 
15.38 17.50 17.97 18.08 18.48 14.36 

Tiiad2 Tvapr Psat Tsat Tsump 
BR 2 0g dos elt S2=-10.5 
Ee 

2.1 14&*03 
Qdo Athetao Ahtude AuQdo 


SS7E+O4 11.171  2.4556*02  3.555E*03 


Bulk 011 2 = 2.0 


3 4 5 8 7 3 
[2.23 (2.735 2.40 .0001.315 239. 39 533.93 
[1 ra L3 E! (S 16 
15:53 ^05. 9899815. 500m os a. Som 6 


Tliad2 Tvapr Psat Tsat Tsump 


Die 2 31 10 8 2.12  -18.8 
Ec 

2.4320E-03 

Qdo Atheta Ahtube AuUdo 


I29E +04 14.092 4.213E+02 3.944E+03 


sulk 011 24210. 4 


j 4 S 8 Hi - d 
[4.14 14.73. 14.52 113.94. 39739 2933 
11 IZ 13 14 IS 16 


(4.70 16.26 15. IMM. 07 AS. > al 


Tliad2 Tvagr Psat Tsat Taump 


2.14 ken 10.70 2.33 -IQ.s 
Ce 

3.972€+03 

ddp Athetao  Ahtube Ruddo 


2.(832E*03 3. 3616+04) 15.302. 12. 755E*02 5. 720E+03 
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Data Jet Number 2 13 3uik 011 X = 2.2 


TIME: 1:35 51553 
TC No: ! 2 S 4 5 5 7 3 
Acne PISA MIGA Ss Ss 7 02 19.37 16.17 -39.99 -39.99 
TÇ No: 3 19 A |: 15 14 IS 18 
Temo : 239.239 15.73 4,44 15.37 TIR 165.93 IS -3 (13.23 
Twa ATwa Da Tliad2 Tvapr Psat Tsat Tsuma 
LEE St eae tA P 2.18 10.37 2.21 -19.9 
Ur le Co 
rA 1.84 JocSdE t3 
Thetab Htuge Qao Athetap  àhtube AuQdo 
14.230 1.182€-03 3.333E-04 13.125 2.353E*02 3.755E-03 


Data Set Numoer = |4 Sulk Oil 1» @.2 


TIME: 13:32:48 

TC No: 1 2 3 4 3 5 7 3 
Tema : 17.77 18.27 17.33 18.48 18.53 17.70 -39.39 -39.99 
TC No: 3 19 11 p ES 14 IS 16 


EROS SS RE ISSO SS SPAS. 2.38 


Tua ATwa Tliad Tliqu2 Tvapgr Psat Tsat Tsumo 
17.54 15.03 2.23 2.35 ZI 282 2.2? 10.9 


Ur ir ge 
3.23 1.26 1.207E*03 
Thetao Htube Ado Athetab àhtube AuQdo 


aes 13. 15559051773. 584E *04 —12.7504 2.3/7 9907  3.7928*03 


Data Set Number = 15 Bulk O11 % = 39.3 


TIME: 1327354720 

Te NO: I 2 5 4 S 3 7 - 8 
Temo : 16.36 17.52 17.20 17.74 17.99 16.30 -99.39 -39.99 
TC No: 9 18 l! 12 13 14 IS 18 


Temo 99.393 15:59 f4.13 15.58 15.539 165.54 15:18 12.34 


Twa ATuwa Tirad Tlrad2 Tvapr Psat Tsat Tsumo 


RS ISS Ss 2.40 2.54 2520. Views 20 Sg Ed t 
Ur le Ce 
3523 0.75 A 
Thetab Htube Qdo Athetab  Ahtuae Ou dp 


14.363 $.S72&+03 65.558E+04 12.362 2.387E*02 3.782E+03 


| 


Data 3et 
TIME: 


WE Nos 
Tamp : 
(C No: 


Numoer = 16 


13:35:45 
I - 


IS JEcas 
a 12 


Temp : -39.39 15.49 


Tua 
15.30 
Vr 
5.38 
Thetab 
13.303 


Data Set 
TIME: 


TC No: 
Temp : 
TC No: 


ATwa Tliad 
14.36 Zea 
De 
L 
dtuoe 
3.431E-*03 


4 


Number = 17 


13:36:51 
1 


13.76 14.46 
3 1Q 


Temp : -99.99 15.21 


Twa 
14.92 
Ve 
SM 
Thetab 
11.348 


Qata Set 
TIME: 


TC No: 
Temp : 
TC No: 


ATua iiio 
14.30 2:28 
Dc 
Is 
Htube 
2.271E+r03 


= 
> 


Number = 18 


13:38:09 


l 2 
117968 12.54 
3 19 


Ténp 3:39.99 209.239 


Tua 

1258 
Ve 
4.26 
Thetab 
10.048 


ATwa Tligd 
14.30 2528 
Ir 
. 87 
Htube 
I .SS8E+03 


l 


Bulk 011 à = 9.2 


3 4 3 3 7 3 
15.38 15.59 [6.32 15.53 -939.39 -39.39 
11 fiz LES 14 1S 16 
14.91 15.51 Ee 

Tliad2 Tvapr Psat Tsat Tsumo 
as 2.2S 19.34 Z.B -19.3 
EE 

3.300E+03 
Qda Athetab Ahtube AuQda 


.334E*04 12.309 3.039E+02 "943 /9E-03 


Bulk Oil X = 0.2 


3 4 5 5 7 3 
14.19. 14.31 . 14.21 015.09 7397390533933 
11 E e 14 is 16 
13.75 135.17 ` IS. IQ 16.15 PiS 0 ame 


Tliqd2 Tvapr Psat Tsat Tsump 
2 d 2.08 10.33 ool ee Oe 
E 

3.910€+03 
Qda Athetab Ahtube AuQao 


sS9QE+OR 12.427 95, 03Sc -0c Sc 34 TECOS 


Bulk Oil X = 2.2 


3 4 S 5 T 
12.42 12.98 12.20 11.69 -39.39 -39.99 
E 12 13 14 IS 16 
13.82 15.29 15.25 16.45 15.37 12.5] 


Tliqd2 Tvapr Psat Tsat Tsumo 
224 2.96 10.32 2.18 -10.1 
Cc 

2.437E+03 
Qdo Athetab Ahtube duQda 


.363E+94 12.638 3.J038+92 3.797E+93 


1952 


Qata 3e 
TIME: 


DG 
Tamp 
TC No: 
Tamo 


No: 


Tua 
3.30 
Ur 
2 72 
Thetab 
ao! 


Qata Se 
TIME: 


TC No: 
Tema 
TC No: 
Temp 


Tua 
Tu] 
Ur 
2.94 
Thetao 
4.3560 


Qata Se 
TIME: 


TC No: 
Temp 

TC No: 
Temp 


Tua 
Su 
Ur 
1.74 
Thetab 
3.336 


t Number = 19 dulk O:!1 < = 2.2 
Ross 
| 2 3 i 3 3 7 3 
3.20 Hee 3.45 3.30 3.35 ssa 
3 19 11 Kä ES 14 1S 15 
-99,39 15.553 14.21 15.148 15.565 16.58 16.32 13.16 
ATwa Pliqd Ttiqul Tvapr Psat Tsat Tsumo 
[3113 ong qua S O IO 2024. —1951 
Tr ES 
«53 1.809E*23 
Htube Qdo Athetao Ahtude AuQdo 
Jg.2496-02 35.712E-*03 12.390 2.368E-02 3.096€+*03 
t Number = 20 Bulk Oil X = 0.9 
13:40:18 
I 2 3 4 S 8 T g 
Z. 12 Tub 7.03 7.10 eg 7.23 -99.39 -99.39 
g 10 11 I I 14 1S 16 
-99.39 15.55 14.29 15.45 15.38 16.07 19.31 13.21 
ATwa ligg Tliad2 Tvaor Psat Tsat Tsump 
SS ENC 2437 uro E091 2.25 -10.1 
Ir Ce 
us 3.373E*02 
Htuba Qdo Athetab Ahtuoe Audo 
7.498E*02 3.717€+93 12.300 2.3I22£r02 3.753E-925 
t Number w 21 Bulk Oil % = 9.0 
13:41:00 
I 2 3 4 S 8 7 E 
8.17 6. 6.34 8.32 5.77 §.33 -39.99 -99.99 
3 19 11 12 13 E! IS 16 
: -99.39 15.78 14.52 15.08 15.39 18.76 15.50 13.48 
ATua Tita Tiiqd2 Tvaer Psat Tsat Tsumo 
IS.34 ERES 2.38 2.25 19.94 2,29 10.4 
Ir Ce 
2o 7.367E+02 
Htube Qdp Athetab Ahtube AuQdp 
5.991€+02 2.577E+03 13.253 2.368E+92 3.744E+03 
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Qata 3e 
TIME: 


TE Nos 
Tamp 
TC No: 
Tamo 


Tua 
d 
Ur 
l; 353 
Thetab 
2.053 


Data Se 
TIME: 


TC No: 
Temo 
TC No: 
Temp 


Tua 
3.58 
Mr 
99 
Thetao 
F S 


Data Se 
TIME: 


TC No: 
Temp 

TC No: 
Temp : 


Tua 
2238 
Ur 
.70 
Thetab 
.374 


NOTE: 


t Number = 22. Bulk Uni O 
13:542 608 
| E 3 4 3 3 7 2 
4.38 3,49 1.39 An 357 4,33 4.31 -39.39 -39.39 
g 10 l l kä l IS E 16 
=39.39 115.38 14.70. 15.95 15.97 15:26 P iS SS MISS 
ATwa Fl qü Tliqd2 Tvapr Psat Tsat Tsumo 
19,3] 2.45 2.231 2993. dd TARAI] 02:9 
Ir Ge 
wen o cres 
atuoe don Athetan  Antube ^uQdo 
7.040E-02 1.572€+93 13.995 2.858€+02 3.7S5E+03 
t Number = 23 Suik Oil %4 =17.3 
13:35:07 
| E 3 4 S 5 d 3 
3.33 3.393 5.52 3.47 3.33 Ses. cessi sms 
g 10 B iz ES 14 1S 18 
-99.39 15.38 14.54 15.88 15.98 18.39 15.54 13.42 
ATwa Tligd Tlrad2 Tvapr Psat Tsat Tsumn 
15.42 2-3 pM 2.19 10.92 2.27 -10.2 
e ES 
20 1.825c-02 
Htube Qdo Athetao Ahtube Ruldo 
3.329E*02 3.597E-02 13.153 2.808E-202 3.395E+-05 
t Number => 24 Bulk Oil % = 2.2 
1555004 732 
| à 3 4 3 3 7 _ 3 
2.36 2.34 2295 E do s 2.99 -99.39 -99739 
g 10 11 Ez i 14 1S 18 
-39.39 15.99 14.63 15.37 16.08 17.01 15.38 13.30 
ATwa Tirad Tliqd2 Tvapr Psat Tsart Tsumo 
15.52 2.48 2.48 2.34 11.95 2.40 -10.9 
Ir Ge 
14 5.344E+91 
Htube Qdo Athetao Ahtube Ou do 
7.392E-02 4.241€+02 13.128 2.947€+92 3.735€+03 


24 data runs were stored 


in file 00907SCZ 
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APPENDIX D. SAMPLE CALCULATIONS? 


Data file DO921ASC4 was used for sample calculations. The 
saturation temperature was 2.2 °C with a heat flux of 59,130 


W/m? on the smooth upper test tube. 


A. SMOOTH TEST TUBE DIMENSIONS 


Do = 0.015875 m 
Di = 0.0127 m 
Die 0 0/2446 m 


ih = 0.2032 m 


Lu = 0.0762 m 


B. MEASURED PARAMETERS 


EE = VroZs —e 236 V)(25) =-194.0 V 


EL— TIr*1.9152 = (1.66 A) (1.9182) = 3.18 A 
Tl = 19.44 °C TNI 2 0599. .^C 
E C m= 19763 C 
T5 = 19.46 °C Ge ISTC 


T7 and T8 not used, defective 


^?This analysis is essentially identical to the analysis 
done by Perry [Ref. 2] and therefore his procedure is largely 
reproduced. 


TSS 


P el 


kc = 344 W/m'K 


Cc. OUTER WALL TEMPERATURE OF THE UPPER TEST TUBE 
1. Surface Area of Smooth Tube 
P = x*Do = x(0.015875 m) = 0.04987 m 


Ac = xr (Do*-Di*)/4 = xw((0.015875 m) "0.01 5 PR 


Ac o ou n 
2. Heater Power 

Oh = Vel = (194.0 V) (3.18 A) =o reo EN 
3. Average Inside Tube Temperature 


Tavg = E," Tn/m = E,* Tn/6 = (T1+T2+T3+T4+T5+T6)/6 


Tavg (19.44+20.03+19.75+19.63+19. 10955000 


Foo 
4. Outside Wall Tube Temperature 
Two is calculated using Fourier's Conduction Law, 


assuming uniform radial conduction and using a known Tavg: 


Two. = Tavg = (Qh(mn (Do/ZDi))/ on ke 

Two = 19.62-(616.92(1n(0.015875/0.01245))/ 
(20S i, 

Iwo ee 


5. Tube Wall Superheat 


0 = Two = Tsat = ]9- 25F- "0 O 0 s Sl sO Rn G 


SiG 


D. R-124 PROPERTIES AT FILM TEMPERATURE 


The 


thermophysical properties of R-124 were computed at 


eee itilm= temperature using the REFPROP equations provided by 


Bertsch 


1. 


[Ref. 1]: 

Film Temperature 

O ES aa IS Oo + 2:20 72 = 10.74 °C 

Dynamic Viscosity 

nuuc 15953 EC EE 
(2207 092105 IDE! 

UI As os O OA) EE 
CEET EE EE EK kat 

m 103421. 0 Ne sy im 

Density 

pd SILO T£=-(6.72617210 i PIS52xI0 ?)Tf^ 

puc 14347 8= (2. eol9) 1074-16. 7267x107) (10774) *- 
EE 

p ^ 1403.2 kg/m” 

Kinematic Viscosity 

u u p = 0422210 / 1403°2.— 2217x10 m/s 

Thermal Conductivity 

F f o> 10 mam5436x10 EE us OCT EE 
5165 5710 Tf" 

k 2 7.5191x107- (3.5436x10*) (10.74) - (1.9545x10"') 
auod E M4 


E3710: W/mK 


ES 


Specific Heat 

Cp — 1.05421(2-1405x10 MESS 
(6.4721x10^)-Tf^-(1.4324x10^) Tf't(4.136x TON DE 

Cp = (1.0542*(2.1405x10) (10.74) * (1:0709x10- 2) (MOSS 
(6.4721x10°®) (10.74) 112492 4: Om elon a 5 
(4.136210 !') (10: 74) NDS 

Cp = 107373 Say Kook 


Thermal Diffusivity 


x k/ (p Cp) = 7137x100: 22 0D 05 2L OIN 


q = 4.717x10^? m/s 

Volumetric Thermal Expansion Coefficient 

B = -Ap/(p-AT£) = - (1402.89-1403.50)/(1403.2:0.2) 
p 


Prandtl Number 


EE 


EECHER 


E. HEAT FLUX THROUGH NON-BOILING TUBE LENGTHS 


Using the Churchill-Chu correlation [Ref. 9] for a smooth 


cylinder in natural convection, the average natural convection 


heat transfer coefficient of one non-boiling tube lenguas 


h 


where 


k/Do{0.6+0.387([F/ (1+ (0.559/Pr) 716) 927] } 


T [[g:B:DO 0 canh (m:bu Jo o r: um) S 


nep ras |” 


m 


Mics pregram DRPJY computes the value of h by assuming a 
starting value of 190 Ww/m^-K and iterating to the final value. 


The resulting values of h and m are: 


h = 281.93 W/m*-K 
m = 23.94 (1/m) 
therefore, 


Of = (h:p:kc:Ac)*?:8 :tanh(m:Lu) = 9.52 W 


F. HEAT FLUX THROUGH BOILING TUBE LENGTH 


Su COn = 20Of — 616702 HL) 597156 W 


J 01 o> 10 mé 


Ab TM bosch = x (0.015875 2032) 


q" = Q/Ab = 597.88/1.013x102 = 59,021 W/m2 
h = q"/8 = 59021/17.08 = 3456 W/m-K 


The following results were calculated by the data 


acquisition and reduction program DRPJY: 
q" = 59,130 W/m? 


h 3463 W/m^-K 


ies SC 


CD 
ll 


E59 


APPENDIX E. UNCERTAINTY ANALYSIS* 


The same data run used for the sample calculations 
(Appendix D) was used for the uncertainty analysis, thus all 
the measured and most of the calculated results still apply. 
A lower heat flux of 1561 W/m’ was also analyzed from the 
same data run to show the difference between the two heat 
fluxes. All uncertainties are presented as a percentagens 
the calculated parameter, i.e. a relative uncertainty. The 


uncertainty analysis method suggested by Kline and McClintock 














(Ref. 29] was used. For example: If R = R(x1,X;,X,, - Ia s se 
then: 
3 Mn ; ns, Cl 
SR EE 5X, 
5x, 5x, Š x; 
where: 
ôR = uncertainty of the desired dependent variable 
X, = measured variable 


dx, = uncertainty in the measured variable 


‘This analysis is essentially identical to the analysis 
done by Perry [Ref. 2] and therefore his procedure is largely 


reproduced. 
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A. HEAT TRANSFER RATE UNCERTAINTY 


T. 


1.66 amps Ir 050/75. amps 


Vr T.16 VOLTS QV T OMS TO e 


so we have: 





ooh . E A t ib 


VE Ir 


B. SURFACE AREA UNCERTAINTY 
The following dimensions were taken from the 


manufacturer's data sheet, and from actual measurements: 


Do 


15:399 (mm) óDo - 0.1 (mm) 


L Zio mm) ôL = 0.1 (mm) 


Ab = Do x L 


ent 


C. WALL SUPERHEAT UNCERTAINTY 


The following temperatures were measured: 


de 
T2 
T3 
T4 
T3 


T6 


M 
240 
1 
19: 
R 
r9: 


44 
03 
13 
63 
46 
38 


ue 
C 
E 
uc 
E 


E 


Tavg = 19.62 


Standard Deviation = 0.245 °C 


Tsat = 2.20 °C and ôTsat = 0.01 °C 
Tw T - Qh | ———— 
ES avg ms SE 
Note: The second term on the right hand side of the above 


equation is the Fourier heat transfer conduction term, and can 


be neglected for the uncertainty analysis since its value is 


less than the standard deviation, so: 


TWO 


AT 


Tavg 


— 
——— 


Two - Tsat 


Dor 02e 
= 19,62 Co 2 20 Cas 1o oe 


1152 


therefore, 


AT 





AAT [|6ATwoy , / ó Tsati2]3 
(ar?) (a) 





BAT [0.2452 ,| 0.01 V] . 
AT vci "E uic 
D. HEAT FLUX UNCERTAINTY 
q"- uno xS 
AD 
where: Qh =VxI = (194 V) (3.18 A) = 616.92 W 
and 60h = (616.92) (0.0164) = 10.12 W 
q" = 59130 W/m^ (measured) 
Ab = 0.010137 m 
Qn - 2.0f = (q)(Ab) = (59130) (0.010137) = 599.40 W 


solving for Of: 


Or j le 02 — 599940)/2 89.76 W 


assuming the same proportion uncertainty for Qf, we have: 


óQf = (8.76) (0.0164) = 0.144 W 


and for the uncertainty: 


tol 





a” ` [es 262) (ao) + EY | 








bq” 105121 g 2 x 0.144) | a E 
S “OO A. 0.0063 = "q" 
ES E | 599.4 | ( ) | 80% 


HEAT TRANSFER COEFFICIENT UNCERTAINTY 








E. 
h = q"/AT 
therefore, 
2 1 
Sd Poa! Ron mao 
h g4 | c 








óh _ 2 nz 
- (0.018) «(0.0141)]? = 2.29% 


B 
The above series of calculations was also conducted with a 
lower heat flux for the same data run and with a high and low 
heat flux from a TURBO-B data run (D1022FTBC4) with the 


results tabulated in Table 5. 
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TABLE 5. RESULTS OF UNCERTAINTY ANALYSIS 


=== Tum. [mmo Dunn [unn] 
[5on/oh (9) |8.49.— — |1.64 — Lin fico —— | 
bAT/AT (3) |3.55 1.41 26762 12.44 

bq"/q" (°) 
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